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INTRODUCTION 



One of the significant developments of modern times 
has been the rapid growth in the size and complexity of 
public and private organizations. Decisions made by mili- 
tary and industrial leaders affect large amounts of capital 
and numerous people. Errors can be extremely expensive and 
an incorrect decision may require years to rectify. Further, 
decisions must be made quickly to match the increasing speed 
and change of modern human endeavors. Failure to act rapidly 
may provide a competitor with a marked advantage. 

Since World War II, men have turned increasing effort 
toward improving the speed and accuracy of decision making 
by application of analytical and objective methods. The 
techniques for optimizing decisions presented in this thesis, 
represent an important part of those methods 

The optimization techniques explored here have potenti- 
ally a wide application in the field of petroleum management. 
These are: linear programming, game theory, and inventory 
theory. Chapter III includes a brief discussion of queueing 
theory and Chapter IV concerns the optimization of purchas- 
ing policies. 

The treatment of each field endeavors to include a sur- 
vey of the extensive literature covering the technique and 
provides an evaluation of appropriate source material. The 
techniques are presented in some detail with specific refer- 
ences suggested where it is considered that the reader may 
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desire exhaustive treatment of a particular facet. Proofs 
are not offered but are generally available in the refer- 
ences. The application of each technique to various mili- 
tary and industrial petroleum problems is handled in the 
form of assumed problems or as a discussion of applications 
believed appropriate to petroleum problems facing the mili- 
tary establishment. 

Before proceeding, it is considered appropriate to 
describe the general characteristics which are common to 
each of the optimization techniques covered in this thesis. 
The methods discussed are those which seek to optimize the 
effectiveness of operating systems. 

Measure of Effectiveness 

In order to compare the effectiveness of various alter- 
native courses of action, a scale of some type, usually nu- 
merical, must be selected as a measure of effectiveness. 

This may be in units of dollars, hours, successes per hundred 
attempts, barrels per day, or it may be dimensionless, such 
as the probability of destruction of a target. The effec- 
tiveness may depend on a number of variables, some of which 
are beyond the control of the decision maker. All factors 
must be accounted for in a model which serves to relate the 
various alternatives to the effectiveness of the system. 

Model Formulation 

Basic to the solution of a problem involving an operat- 
ing system, is the development of a model. Models can be 
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highly mathematical - and either deterministic or probabalis- 
tic - or they can be qualitative and non-mathematical. The 
primary function of a model is to provide a means of bring- 
ing together the facts and data bearing on the problem. As 
a substitute for the real system, the model can be manipu- 
lated to test various alternatives in a way which an actual 
system can rarely be handled. The degree of correlation 
between the actual system and its model is a partial measure 
of the limitations which must be recognized in the solution. 

Sub-Optimization 

It is possible to optimize the cost of operation of a 
catalytic cracking unit but to do so at the expense of the 
overall long term profits of the refinery. The selection 
of the proper level at which to optimize operations, is an 
important decision. Optimization at too low a level leads 
to sub -optimization. At the same time, the problem at a 
higher level may become so complex that an adequate solu- 
tion becomes extremely difficult if not impossible. 

The optimization techniques presented in this thesis 
include an important segment of the analytical and objec- 
tive methods which are being -increasingly applied by govern- 
ment and industry to achieve better decisions. 



CHAPTER I 



LINEAR PROGRAMMING 
BACKGROUND 

The techniques of linear programming have been utilized 
extensively and with considerable success in many instances. 
It should be noted that such techniques are applicable to 
some, but not all of the problems encountered in programming. 
Among the numerous operating systems, there are many in 
which the fundamental relationships of the variables are 
linear in character. The scientific study of such systems, 
with the aim of optimizing selected factors, has led to the 
development of the theory of linear programming. 

A number of authorities have indicated that linear 
programming is related to and developed from economic theory. 
Historically, VonNeumann, Neisser and others recognized 
during the 1930's that the simple version of the Walrasian 
general equilibrium could not be adequately treated in terms 
of the number of equations and unknowns. More recent inter- 
pretations of the early work of Frederick Taylor, Smith and 
Walras by economists such as Samuelson, gave impetus to the 
development of linear programming. 

In 19^5* George Stigler developed a solution to the 
diet problem which involved determination of adequate diet 
at minimum cost for 77 foods and 9 nutrients. He determined 
that optimum diet in 1939 consisted of wheat flour, cabbage 
and dried Navy beans and cost $39.93 Tor the year 1 . 
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Organized research in the area of linear programming 
was conducted in the U. S. Air Force Project SCOOP (Scienti- 
fic Computation of Optimum Programs). The outstanding con- 
tribution of this project was made by George B. Dantzig in 
19^7 when he set forth the mathematical statement of the 
general linear programming problem and developed a means 

t 

of solution known as the simplex method. Prior to this 

time, many problems had been unsolvable even though recog- 

p 

nized as linear programming problems . 

Other contributions were made by Hitchcock (19^1) and 
Koopmans (19^7) who independently developed the transpor- 
tation problem and by Charnes who extended the application 
of the simplex method. 

The difficulty of solution of long linear programming 
problems by manual calculation has led to the use of elec- 
tronic computers which add nothing to the theory but pro- 
vide tremendous calculating speed. The first successful 
use of a computer in the solution of a linear programming 
problem was on the Bureau of Standards SEAC in 1952^. 

Most electronic computers have now been programed to solve 
the simplex algorithm and linear programming has become a 
powerful tool in the optimization of operating systems. 

EVALUATION OF BIBLIOGRAPHY 

A very real problem exists in selecting sources of in- 
formation for preparation of a paper on linear programming 
methods and applications. While the material available on 
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the subject has sprung forth only in the years since World 
War II, the volume of articles, books and studies in this- 
area has become a veritable flood. It has not been possible 
to review more than a portion of the books and articles in- 
volving linear programming, but a conscientious effort has 
been made to select’ those which appear to be applicable in 

the light of my subject. The 865 page annotated bibliogra- 

4 

phy. Operations Research by Batchelor , has been of consider- 
able assistance and is recommended to those seeking sources 
on methods and applications. 

The basic work in linear programming done by Dantzig, 
Charnes, Koopmans, and Kuhn has led the author to utilize 
works by these men in the theoretical areas where possible. 
However, the developments of these men are effectively pre- 
sented in more readable form in the basic text covering 

e; 

linear programming by Gass- / . For those interested essenti- 
ally in learning and applying the techniques of linear pro- 
gramming, the book by Gass is recommended as well as the 

r 

sections on linear programming by Saaty and Churchman, 

Ackoff and Arnoff^. 

For applications to military and industrial petroleum 

8 Q 

problems, the work done by Alan S. Manne * , and Gifford 
H. Symonds^ is excellent in the areas of refinery operation. 
For applications involving tanker scheduling and other as- 
pects of the transportation problem, the work done by 
Dantzig and Fulkerson^ as well as Flood"^, provides an ex- 
cellent basis for further study. 
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For the reader interested in additional research in 

the area of linear programming, it is suggested that a 

review be made of the bibliography on linear programming 

ix 

by Vera Riley and Saul Gass 

Mathematical Methods of Linear Programming 

Programming is the planning of an operation for the 
purpose of optimizing the result. The objective may be to 
minimize costs or to maximize returns. It may also be to 
allocate efficiently limited or constrained resources to 
meet specific objectives. The number of solutions may be 
infinite, but the best or optimum solution must satisfy 
both the conditions or constraints of the problem as well 
as the objective. 

When linear (straight line) constraints are assumed 
together with a linear objective, optimization requires the 
solution of a linear programming problem. Assumptions of 
linearity are frequently appropriate, however, it may be 
found that the objective function is non-linear and that 
use of such a function improves the result. It is considered 
that non-linear programming, is outside the scope of this 
thesis and no treatment of the subject is included. Exten- 
sive information on this subject may be found in 14 and 15. 

The General Linear Programming Problem 

The mathematical model of the general linear programming 
problem consists of a set of simultaneous linear inequations 
which represent the conditions or constraints of the problem 
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and a linear function which describes the objective of the 
problem. Mathematically a linear programming problem may 
be formulated by letting a^, b^, and c^ be sets of known 
coefficients with (i = 1 , . .., m; j =1, . .., n) and also 
letting Xj(j =1, ...,n) be a set of unknown variables. 

In the general linear programming problem, we seek 
solution sets which may be regarded as vectors, X = (x^, x g , 
..., x n ), which satisfy the linear inequalities previously 
referred to as constraints, and at the same time maximize 
the linear objective function. 

The inequalities in the general case are of the form: 

n 

a^j Xj ^ b^ (i = 1, . . . , m) (l— l) 

j=l 

The foregoing is a condensed notational form which can 
be more easily understood in expanded form as shown here: 

a ll X 1 + a 12 x 2 + • • • + a i n x n S b i 

a 21 X 1 + a 22 x 2 + • • • + a 2n x n £ b 2 



a m ,3, + a 0 x 0 + . . . + a mv x £ b m 
m, l l d mn n — m 

It is also established that the values of the variables 

are positive and the solution set is thus subject to the 

constraint: 



x 3 s 0 



(j — 1, 



n) 



( 1 - 2 ) 
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The linear objective function which is to be maximized 
is of the general form: 



The objective function can be expanded in the more 
easily understood form: 

Z = C^x^ + ^ 2 x 2 + ♦ • • + C j x j 

Before attempting to illustrate a solution, let us 
examine the problem structure briefly. Equation (1-5) de- 
fines the objective function, Z, which is linear in each 
value of Xj . It is a function of the vector X = (x^, x^, 
..., x ) and may be expressed as: 



This function is defined for all values of x with finite 
components, but must also satisfy the restrictions of equa- 
tions (l-]) and (1-2). 

Equation (1-2) is usually called the non-negativity 
restriction and requires that only positive values of x be 
considered. Equation (l-l) requires that feasible solutions 
satisfy m linear inequalities generally termed "constraints". 

The linear programming problem may be stated in an 
equivalent form where a linear function is to be minimized 
and the constraints are ^ instead of Since this problem 
is obtained from the general problem described by multiply- 
ing the inequalities by -1 and maximizing -Z, this case can 
be covered by limiting this discussion to the maximization 
problem. 




(1-3) 



j=l 



Z = f (x) 



(1-4) 
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Methods of Solution 

The primary methods of solution of the general linear 
programming problem are the geometric method, the relaxation 
method, the method of double description and the simplex 
method. Consideration of the geometric method will serve 
to illustrate and clarify the problem and is therefore in- 
cluded at this point. A discussion of the simplex method, 
which is the most frequently used technique, will also be 
given. 

Geometric Method 

Geometrically, a solution to the linear programming 
problem is a point of the convex set defined by the con- 
straints which also maximizes (minimizes) the objective 
function. 

To illustrate, let us consider a simple production 
problem of a manufacturer of two types of lubricating oil. 
Both items require the use of common equipment and labor. 

The objective of the company is to maximize profits on the 
two products. It is assumed that a fixed market price is 
established and that there are six fixed resources to be 
allocated to the production of the two items. These are 
the hours of skilled and unskilled labor available, the 
number of machine hours available on two blending machines 
through which both products must pass and two items of 
blending stock required as raw materials. 

Let x^ and Xg be the respective quantities of the two 
lubricants to be produced. Assuming 5 and 8 to be the 
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respective unit profits, our goal is to maximize the objec- 
tive function 

Z = 5x 1 + 8x 2 



The constraints are expressed as: 



(1) x 2 £ 5 

( 2 ) Xj + 2 x 2 £7 

(3 ) 2x 1 + x 2 < 8 

(4) 3x 1 + x 2 < 13 

(5) 4 X;l + 3x 2 < 40 

(61 3x x + 2 x 2 < 20 



Skilled labor 
Unskilled labor 
Machine 1 
Machine 2 
Blending Stock A 
Blending Stock B 



and non-negativity restrictions are 



*i s 0 

x 2 

Because of the non-negativity restrictions we are con- 
cerned only with the first quadrant and the planes of the 
constraints appear as lines in two dimensions. See Figure 1. 
The solution set which satisfies the conditions of the con- 
straints is thus a polygon as shown and is both bounded and 
convex. We are thus assured of at least one optimal solu- 
tion and there are alternate optimal programs only if the 
objective function is parallel to one of the sides of the 
polygon. It will be noted that the only constraints which 
affect the solution set are those of skilled and unskilled 
labor and machine hours on the first machine only. All 
other constraints are redundant and of no effect. 
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Figure 1 

Set of Solutions 

Figure 2 illustrates the geometric method by consider- 
ing a family of planes for different values of the objec- 
tive function. It will be noted that the maximum Z, being 
the greatest distance from the origin, is a line through Pg 
with Z = J1 and a slope of - 8 / 5 . The optimum solution thus 
calls for production of 2 units of Xg and 5 units of x^. An 
analysis of the solution shows that even the hours of skilled 
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Figure 2 

The Optimum Solution 

labor available is not limiting, and further savings could 
be made by adjusting inventories and resources to the end 
that all constraints will pass through the point P 0 , the 
optimum solution. 

Simplex Method 

The simplex method as developed by Dantzig provided the 
first practicable and simple system for solving linear 
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programming problems with large numbers of variables. As 
originally proposed the method did not adequately handle 
certain degenerate cases; however, work by Charnes resolved 
these problems and the simplex method was made workable with- 
out the restriction of non-degeneracy. The most complete 

* 

coverage of the simplex method is set forth by Charnes . 

This treatment is exhaustive. For the reader interested ' 
in simply learning to apply the technique, the presentations 
given by Gass' 1 , Churchman -1 , and Sasieni y are recommended. 

The simplex method is an iterative procedure and each 
iteration consists of translating the plane of the objective 
function parallel to itself for successive values of Z and 
evaluating its distance from the origin at the vertex of the 
convex set. This can be more adequately explained by re- 
ferring to figure 2. Although the figure is in only two 
dimensions, the convex set can be visualized as describing 
a polyhedron formed by the planes of the constraints. One 
seeks the point of the polyhedron (generally on the boundary 
and usually at a corner or vertex) through which the Z plane 
may pass and which plane corresponds to a maximum value of 
Z. In the simplex process planes of Z are passed succes- 
sively through the corners or extreme points of the polyhed- 
ron, advancing away from the origin with increasing values 
of Z, until a maximum (optimal) solution is obtained. 

The simplex process also utilizes a criterion which 
eliminates some of the points as possible solutions. The 
number of iterations required varies but is frequently about 
twice the number of inequalities. 
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The first step in the simplex method is to convert the 
system of m linear inequalities to a system of m linear 
equations. This is accomplished by. the insertion of "slack 
variables", x x^g* ••**> x n+m' one l‘ or each of the in- 
equalities which then becomes: 



a ll x l + 


a 12 X 2 + 


. . . + 


a ln x n + x r>H 


* 


“ b l 


a 21 x l + 


^22^2 “t** 


... + 


a 2n x n + + x n+2 




" b 2 



(1-5) 



a.x, + a ~x 0 + ... + a^x^ , . 

ml 1 m2 2 mn n + x . = b 

n+m m 

The use of the slack variables above is simply a device 

used in arriving at an optimal solution. The use of these 

slack or idle factors will be more readily understood in 

the numerical example to be given later in this chapter. 

Having converted the inequalities to equalities (1-5), 

the problem can be restated as one of finding among all sets 

of points (x^, x^, ...., x n+m ), that point which satisfies 

the constraints (1-5) with non-negative values of x. and 

J 

also maximizes the objective function: 

Z = C^x^ + C 2 x 2 + ... + C n x n (1-6) 

The basic theorem behind the simplex technique states, 
in effect, that a set (x-,, x,., ...., x ,) which maximizes 
the objective function and satisfies the constraints, must 
necessarily have n elements which are zero. If the coordi- 
nates were known, we could substitute their values in the 
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constraints and would then have m equations and m unknowns 
which could be solved algebraically for the optimal solu- 
tion. Without this knov/ledge, it is necessary to proceed 
step by step to a solution. 

The procedure may be explained as follows: 

Choose n of the n + m variables and give them a value of 
zero. Utilizing these values, the constraints may then be 
solved for the remaining m variables. Solutions with nega- 
tive values for some variables are discarded as violating 
the non-negativity restrictions and another set of n vari- 
ables is selected for trial as zero variables. When the 
values of all coordinates or variables of the solution set 
are non-negative, we have a feasible solution. This is the 
starting point of the iterative procedure. 

The simplex method also permits a test of the feasible 
solution for optimality, and, if not optimal, the method 
allows one to proceed directly to a new feasible solution 
with an improved Z. Ultimately the procedure leads to an 
optimum solution. 

Application of the simplex technique involves extensive 
use of matrix representation and manipulation of the rows 
in a series of steps until a solution can be read from it. 

To illustrate the manner in which this is done, a numerical 
example will be presented here. 

Numerical Example of Simplex Method 

In this example, the goal is to maximize the objective 



function: 
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Z = 3x 1 + 5*2 + (1-7) 

subject to the linear constraints: 

+ ;5x 2 < 8 

2x 2 + 5x^ < 10 (1-8) 

3x 1 + 2 x 2 + 4x^ < 15 



Step One. Use of Slack Variables 

By introduction of "slack variables", which essentially 
are additional variables, the inequalities become equalities 
The slack variables are non-negative and take the form, x^, 
x^, and Xg in this example. The inequalities (l-8) thus 
become the following equations: 



2x 1 + 5x 2 + x^ =8 

2Xg + 5x^ + x 5 =10 

3x 1 + 2 x 2 + + Xg = 15 



(1-9) 



The first goal of the simplex method is to obtain a 
feasible solution associated with an extreme point of the 
solution set. Such a feasible solution is available after 
the first "tableau" or matrix of coefficients has been 
formed. However;, it will be necessary to continue to exa- 
mine other extreme points in a series of tableaus until a 
maximum solution is obtained. 

Equation (l-9), can now be written in matrix form, mark 
ing the columns of coefficients as P^, ... Pg and using 

P Q for the column of constants. Below the matrix is row 
where the coefficients of the objective function are entered 
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First Tableau 





1,2 



t 

entering 

variable 



8/5 10/2 15/2 

departing 

variable 



Step Two. First Feasible Solution 

In the foregoing matrix, a feasible solution can be 
read off, so long as the non-negativity restriction is not 
violated. This solution has been entered in a row below 
the objective function. 

This first feasible solution was obtained by setting 
the variables x^, Xg, and x^ equal to zero. The solution is 
read off by selecting the non-zero variables to be those 
whose columns contain a single "1", with the rest being zeros. 
Columns P^, P,_, and Pg are thus chosen and for convenience 
are shown in the P 1 column in the box at the left. The co- 
efficients of corresponding variables, x^, in the objective 
function are inserted in the column marked in the same 
box, again for convenience. With all variables at zero in 
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the constraining functions, except the slack variables, the 
feasible solution is x^ = = x^ =0 and x^ = 8, x^ = 10, 

xg = 15- This is indicated in the solution row of the first 
tableau. The value of Z associated with this solution can 
be seen to be- zero' (Z =»0) since the variables x^, and 

x^ in the objective function have values of zero for this 
first feasible solution. 

Step Three. Test for Optimality 

This is accomplished by computing an evaluation, A^, 
for each zero variable, x^, in the solution by the formula: 

i=m 

A j=°j-I a ij c i (1 - 10) 

i=l 

At this point in the process, where the set of non 

zero variables is just the set of slack variables, it can 

be seen that A. = C. for each j. This can be explained by 
J J 

pointing out that in computing A. for the first variable in 

J 

the solution of our example, we take the value of which 
is 5 and subtract from this the sum of the products of the 
three values in the column times their corresponding 
values in the column designated P 1 in the matrix. 

Since all values in the column are zero at this 
stage, the products are zero. Thus the summation term is 
zero and: 

= 5 - 0 = 5 
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Similarly A^ = 5 and A^ = 4. The A values are entered in 
the row below the solution in the first tableau. 

Interpreting the Test 

If one or more of the A. are positive, the solution is 

J 

not optimal and we must proceed to the next step. If all 
of the Aj's are zero or negative, the solution is optimal. 
However, if in the latter case, some of the A^.'s are zero 
while others are negative, more than one optimal solution 
exists having the same value for Z. If all A^'s are nega- 
tive and none are zero, a unique optimal solution has been 
found. 

Step Four. Computing the Entering Variable 

To form the next matrix, we must find the entering 
variable, x^, which is that zero variable in the feasible 
solution which is to be non-zero in the next iteration. 

The entering variable must have the following character- 
istics: 

1. It is zero in the present matrix 

2. It has a Aj ^ 0 

3* One of the matrix coefficients, a.., in the column 

J 

x. must be greater than zero 
J 

While several variables in a solution might be chosen, 
the problem solution converges quickly if the entering vari- 
able, x r , is chosen as that one £or which is largest. 

In checking the first solution it can be seen that the 
second variable in the solution, has the necessary prerequisites 
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and has the largest A^. It is thus marked "entering vari- 
able" in the first tableau. 

Step Five. Determining the Departing Variable 

The departing variable, x , is that non-zero variable 

o 

in the feasible solution which is to become zero in the next 
iteration. The departing variable is generally determined 
by the selection of the entering variable and by the non- 
negativity restrictions on the coordinates. It is first 
necessary to compute a quotient: 

V a ir (1 ' 1:L) 

for each non-zero term in the solution. The values of b^ 
are located in the solution while the values of a ir are 
found in the column of the previously selected entering 
variable. For the first quotient we compute: 

b i/ a ir = | 

From among the three quotients calculated, the one in column 
P^, is the smaller and thus determines the departing vari- 
ables. The departing variable is designated in the first 
tableau. 

Step Six. Calculate the New Matrix 

To calculate the new coefficient matrix, it is neces- 
sary to perform several rox* operations on the prior matrix. 
The need for the entering and departing variables now be- 
comes apparent. These two variables locate a position in 
the matrix which is marked with a box. The new matrix is 
to have a 1 in this position. By dividing the elements of 
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the first row by 3, which is the value in the designated 
position, unity is obtained in that position in the new 
matrix shown here. 



2/5 


1 


0 


i/5 


0 


0 


8/5 


0 


2 


5 


0 


1 


0 


10 


3 


2 


4 


0 


0 


1 


15 • ■ 



The new matrix also must contain zeros in the other 
positions of the designated column. This is done by sub- 
tracting appropriate multiples of the new row from the 
other rows of the matrix. 

The second intermediate matrix then becomes: 



2/5 


1 


0 


i/5 


0 


0 


8/5 


-4/5 


0 


5 


-2/5 


i 


1 


14/5 


5/5 


0 


4 


-2/5 


0 


1 

1 


29/5 

1 



The new tableau can now be written. 



Second Tableau 

°i p i P 1 p 2 p 3 p 4 p 5 P 6 p 0 
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0 


i/5 


0 


0 


8/5 


-4/5 


0 


5 


-2/5 


1 


0 


14/5 


5/5 


0 


4 


-2/5 


0 


1 


29/5 


i 5 


5 


4 


0 


0 


0 


0 



Solution 0 8/5 0 0 14/5 29/5 (Z=40/5) 



a 3 -i/5 

V a ir 



4 -5/5 

r - 14/15 29/12 

T 4 

entering departing 

variable variable 



00 
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Here it will be noted that the second feasible solution 
is again read off by selecting the non-zero variables to be 
those whose columns contain a single "1" with the rest zeros. 
This marks the beginning of another iteration. By repeating 
steps 1 through 6 , an optimum solution will be reached. In 
this problem, a total of four tableaus are required in order 
to reach an optimum solution. The third and fourth are 
obtained in the same manner and are shown below. 



Third Tableau 



C i p i p l p 2 P 3 P 4 ' p 5 p 6 P 0 



Solution 0 8/5 14/15 0 0 89/15 (Z=256/l5) 

Aj 11/15 - - -17/15 -4/5 
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0 1/5 
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-4/15 0 


1 -2/15 


1/5 
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14/15 


41/15 o 


0 - 2/15 


-4/5 


1 


89/15 



3 5 4 0 0 0 



b i/ a ir 



t 

Entering 

variable 



89/41 

I 

Departing 

variable 



4 



-7/2 
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Fourth Tableau 



c i 


P i 




p i 


P 2 


P 5 


p 4 


P 5 


P 6 


P 0 


5 


P 2 




0 


1 


0 


15/41 


8/41 


-10/41 


50/41 


4 


P 5 




0 


0 


1 


-6/41 


5/41 


4/41 


62/41 


3 


P 1 




1 


1 


0 


-2/41 


-12/41 


15/41 


89/41 




c 


1 


' 3 


5 


4 


0 


0 


0 




Solution 89/41 


50/41 


62/41 


0 


0 


0 (Z =76 5/41) 




A 


i 


- 


- 


- 


-45/41 - 


-24/41 


-11/41 





Since in this final tableau, all A^ are non-positive, 
an optimum solution has been obtained. 

Degeneracy 

The problem of degeneracy exists in the simplex tech- 
nique when the vectors chosen are not linearly independent. 
In such cases, more than n vectors pass through the same 
point which in turn is described by the fact that more than 
m of the original m + n variables are zero. Under these 
conditions the value of the objective function remains un- 
changed by each iteration of the simplex method and the pro- 
cedure continues to cycle. 

Degeneracy problems seldom occur in actual cases and, 

therefore, no attempt will be made to cover the means of 

20 

solving such problems here. Gass has an extensive discus- 

21 

sion of the solution of degeneracy problems and Saaty 
describes a technique for solution of degeneracy cases. 
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Duality 

One other aspect of linear programming should be men- 
tioned before proceeding to the types of problems and appli 
cations possible with linear programming. This is the 
presence of a dual problem which is related to the original 
problem both in formulation and in the existence of a solu- 
tion. 



The pair of linear programming problems shown here are 
duals of each other: 



Find X — (x^> . . . . y x n ) 
such that: 

x i ^ 0 

I a U x j s b i 

J 

^ Cj x j = Z is a maximum 

j 



Find Y = (y^ , y m ) 

such that: 

y ± ^ 0 

I a u y i i c j 

i 

X ^i^i = w is a minimum 
i 



The coefficients a^j, b^ and Cj are constants with 

equal values in both the primal and the dual. Under these 

conditions it is apparent that there is a close relation 

22 

between the two solutions. Kuhn and Tucker indicate the 
following properties are among those characteristic of dual 
programs : 



a. Either both the minimization and maximization have 
optimal solutions or neither does. 
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b. A feasible vector X is optimal, only if there is 
a feasible vector Y such that Z = W. Similarly, 
a feasible vector Y is optimal if there is a 
feasible vector X such that Z = W for these vectors. 

c. If both the primal and dual problems have optimal 
solutions, then max Z = min W. 

These properties indicate that the maximum value of Z 
can be found by determining the minimum value of VJ, or vice 
versa. If W has no minimum value, it then follows that Z 
has no maximum. Thus in any problem in which only min W or 
max Z are to be determined, either the primal or the dual 
may be solved. 

In selecting the problem to be solved, primal or dual, 
the one most easily solved should be used. If geometric 
techniques are used, the problem with the least number of 
variables should be chosen. The reverse is generally true 
for the simplex method which is well suited to the solution 

px 

of the dual. Churchman ^ indicates that the number of itera- 
tions required for the simplex method is from 1 to 1 l/2 
times the number of rows or constraints. Examples of con- 

Oil 

version from primal to dual are given by Gass 

Special Linear Programming Problems 

Of the many forms of linear programming problems, several 
can be grouped by type and special techniques have been 
developed for their solution. Among these are the assign- 
ment and transportation problems which are important enough 
to discuss briefly here. 
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The Assignment Problem 

This group consists of problems where there are n re- 
sources and n uses to which the resources are to be applied. 
Each resource or origin is to be associated with one and 
only one use or destination and we wish to assign the asso- 
ciations in such a way as to maximize (or minimize) the 
summed effectiveness. Problems may involve supplies, person- 
nel, or equipment to be assigned to demands or jobs. Among 
those who have made significant contributions in this area 
Hitchcock, Flood, Kuhn, and Dwyer. The techniques have been 
applied profitably in making optimal assignments of person- 

OR 

nel, equipment, etc. A bibliography by Riley and Gass 
contains references to a number of such applications. As 
a special case of linear programming, the assignment problem 

p 

can be solved by simplified means. In an n matrix, if 

n = 8, there are 8.' or 40,320 combinations to be examined. 

A simplex solution could be used but a technique developed 

by Flood is simpler and can be utilized without use of a 

computer. A detailed example of this method is given in 
26 

Flagle . 

The Transportation Problem 

The transportation problem is a generalization of the 
assignment problem. The matrix of effectiveness is no 
longer necessarily square; however, the problem is essenti- 
ally one of selecting optimum programs for distributing a 
homogeneous product where s^ units are positioned at origin 
0^ and d. units are required at destination D^. The cost 
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of transportation of a single unit from 0 i to is denoted 

by and our objective is to minimize C = ^ x j_ j 

ij 

the total cost of distribution where supply equals demand. 

If supply is less than demand and no preferential destina- 
tions exist, then we can assume the availability of fictici- 
ous supplies which can be distributed at no cost. This pro- 
cedure makes it possible to treat the problem as if supply 
equaled demand. The variable x^ represents the quantity 
of units moved from 0 i to and therefore 

x.j. = s^ and ^ x^. = d^. The condition of supply equals 
i i 

demand is denoted as 



Transportation problems can be solved by the simplex 
method and other techniques but the "transportation technique" 
is the one most frequently used. Several examples of the 
transportation problem are included in most texts on opera- 
tions research. An extensive transportation problem solved 

“ ' 27 

by the simplex method can be found in Saaty 1 . 



APPLICATIONS OF LINEAR PROGRAMMING 
The JP5 Jet Fuel Problem 

A question which is currently a concern of the U. S. Navy 
is one regarding the capacity of the petroleum industry to 
produce JP5 jet fuel. This information is desired, not only 
to determine the availability of JP5 during periods of war 
or peace, but also to consider expanding use of this type of 
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jet fuel and the effect that increased demand may have on 
the price to be paid for JP5. 

No attempt will be made to answer these questions 
directly, but an approach to the problem will be discussed 
as an application of the use of linear programming tech- 
niques to solve military petroleum problems. A survey of 

attempts to answer similar questions, reveals that the work 

28 

done by Dr. Alan S. Manne in 1956 for the Rand Corporation 
in developing a linear programming model of the U.S. petro- 
leum refining industry, can well serve as a basis for solu- 
tion of the questions concerning JP5. The mathematical 
model, constructed by Manne, was one of a series of process 
analysis models aimed at depicting the technological capabi- 
lities of the entire national economy. In general, the pur- 
pose of these studies was to provide quantitative answers 
on production capabilities and also to weigh the effects of 
product substitutions and alternate production processes. 

The general question Manne sought to answer was: With 

the crude oils and refineries available in the U. S. , what 
product-mix possibilities existed between the output of JP4 
jet fuel and the output of other refinery products, and 
what would be the effect of a loss of a portion of refining 
capacity? 

To answer these and other questions in regard to JP5 
jet fuel requires the construction of a linear programming 
model covering both U. S. crude oil production and the opera- 
tion of the nation's refining industry. The model must 
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incorporate the details of the manufacture of JP5 and all 
other major products: 115/1^5 avgas; 100/150 avgas; JP4 

jet fuel; kerosene; premium mogas; regular mogas; diesel; 
Navy special; the several grades of fuel oil; lube oils; 
asphalt; coke; liquefied petroleum gases; and basic aroma- 
tic chemicals. 

The JP4 jet fuel model developed by Manne, did not 
take into consideration the problems of geographic loca- 
tion or the limitations imposed by transportation, new 
facilities or accumulation of inventories. It is under- 
stood that further work in this area has been done by T. A. 
Marschak but information on this work could not be located 
at this writing. Manne indicates that Mahschak has stream- 
lined the technological details of the model discussed here 
and constructed a four region model, (East, Midwest, Gulf 
Coast, and West) which includes limitations on transporta- 
tion. Such a model would probably give an improved repre- 
sentation of production, refining and transportation. Manne 
acknowledges that a model including provisions for location 
and time can be constructed but that the costs of computa- 
tion were prohibitive at that time (195^). 

The model constructed by Manne based calculations on a 
reference date, 1 January 1953 , and assumed that both equip- 
ment capacities and crude oil availability were set at the 
maximum rate estimated for that date. He based all specifi- 
cations (octane number, performance number, boiling range, 
aromatics content, vapor pressure and viscosity) on the 



31 



averages prevailing on that date. He further assumed that 
labor, power, catalysts and sulfuric acid were not limit- 
ing. The model did, however, consider the refinery inputs 
of isopentane, natural gasoline, tetraethyl lead, heat, and 
C2, CJ and C4 gases in addition to crude oil. 

One of the simplifications in the model was considera- 
tion of 25 categories of crude oil only. Based upon data 
compiled by the U. S. Bureau of Mines, a sample of the 25 
leading U. S. crude oils was selected. Because of the vari- 
ations in analysis of individual crude oils, assuming a 
single "composite" crude would not allow for the tendency 
of various types of refinery equipment to utilize the most 
suitable crude. 

It was also assumed that the various types of refining 
processes can be designated by broad categories but not by 
individual patents. No distinction was made between Fluid 
Catalytic Cracking of Universal Oil Products Company and 
the Thermofore Catalytic Cracking of the Houdry Corporation 
since product yields of the two processes are similar. 
Catalytic Cracking and Catalytic Reforming, however, are 
distinct and were treated separately. 

The Linear Programming Model 

The model developed by Manne was a conventional linear 
programming type with numerous constraining inequalities and 
an objective function to be maximized. The model allowed 
numerous possibilities for varying the product-mix. Choices 
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could be made at the three main steps of the refining se- 
quence: atmospheric crude distillation, conversion of 

straight-run products into blending stocks and the blend- 
ing of the end items. In general the constraints intro- 
duced in the model were in six classes: (l) equipment capa- 
city, (2) crude oil availability, (J) refinery gases and 
straight-run streams, (4) converted streams, (5) gasoline 
and jet fuel specifications, and (6) end item requirements. 

The objective function to be maximized represented a 
product mix which, except for JP4 in this instance, was 
proportional to the actual 1955 production of each item. 

It was the purpose of the model, then, to maximize the 
standard product mix, subject to producing specified amounts 
of jet fuel. The requirement for jet fuel can be varied in 
the model to coincide with the expected needs of the mili- 
tary in a given situation. The model, in this way, provides 
a substitution curve between the one group of refinery pro- 
ducts and the jet fuel. Manne referred to these as "trade- 
off" curves. The matrix for the model consisted of 105 rows 
and 310 columns and solutions were obtained by the simplex 
method which was programed for computer use. No attempt 
will be made to reproduce the details of the model or its 
solution here; however, it may be informative to show the 
"trade-off" curves obtained and to discuss something of its 



interpretation. 



Product -mix, exclusive of jet fuel (millions of barrels per calender day) 
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JP-4/jet fuel (millions of barrels per calender day) 

Figure 3 

U. S. jet fuel trade-off curves 
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Figure 4 



Comparison between market prices and shadow prices 






Figure 3 represents the trade-off curve between JP-4 
jet fuel and the standard-product mix - with the various 
assumptions and simplifications given It will be noted 
that, except for-. short segments, the relationship is not 
linear. The curve indicates that jet fuel cannot be sub 
stituted for other products at a constant barrel for barrel 
rate nor at a constant dollar for dollar rate. Examination 
of the curve shows that as the volume of jet fuel increases, 
the amount of other products which must be sacrificed be 
comes increasingly greater. At low volumes of jet fuel, the 
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trade-off curve shows that each barrel of jet fuel can be 
produced at a cost of about 0.6 barrels of standard product 
mix. However, at the extreme of about 5-29 million barrels 
of jet fuel, 25 barrels of product mix is sacrificed for 
each barrel of jet fuel. These figures should not be con- 
sidered as exact amounts. Manne indicated that the curves 
could only be considered- as estimates. In actual circum- 
stances, industry would make adjustments and innovations 
not possible to include in the model. 

The point shown as actual output serves to indicate 
that the model came fairly close to estimating actual out- 
put. However, there are some factors within the model, 
principally the supply of crude oil, which prevent the esti- 
mate from straying too far from the actual. Manne instituted 
an additional check on the model by comparing the shadow 
prices of end items with actual market prices at the time. 

A shadow price, in this model, represents the cost of an 
item measured in terms of the standard product-mix. Given 
the assumptions of the model and neglecting market factors 
(which certainly are present for individual products), 
prices should be proportional and lie on a 45° line between 
the two price scales As ,will be noted in Figure 4, while 
the prices were not perfectly proportional, the correlation 
was quite good and rough estimates of prices may be obtained 
so long as the prices selected for comparison are within a 
reasonable range of actual production volume. Such a pro- 
cedure applied to a model designed to predict production of 
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JP5 by the refining industry, may provide a rough estimate 
of prices although it must be recognized that errors are 
inherent in such a procedure. 

In a second set of calculations made by Manne, it was 
assumed that a reduction had been arbitrarily made in the 
amount of refining equipment available. The trade-off 
curve for the reduced equipment conditions is plotted in 
Figure 3. The linear programming calculations and the 
matrix were unchanged except for the reduced equipment con- 
ditions. Capacities in this instance were reduced to a 
range from 38 per cent to 51 per cent of the earlier levels. 
The purpose was again to maximize the level of the non- jet 
fuel product-mix subject to specific requirements for jet 
fuel. As before, the first increment of JP4 was obtained 
at small cost in terms of other fuels, but it increases 
rapidly to about 17 barrels of product -mix per barrel of 
jet fuel at a production rate of 1.7 million barrels of JP4. 
Prior to the computation, it had been expected that at re- 
duced capacity, it would be possible to select more desir- 
able crudes and a much larger portion of the product -mix 
could be converted into jet fuel. However, there was only 
a relatively small difference in the maximum fraction of 
jet fuel which could be produced from a barrel of crude. 

In the initial calculations this fraction reached a maximum 
of 45 % and with reduced equipment the fraction of crude con- 
verted to jet fuel had a maximum of about 50 %. A more 
striking difference had been expected. 
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It is considered that a similar model can be con- 
structed to provide answers to questions on production of 
JP5* The requirements made in the model by Marschak to in- 
clude the effects of location and transportation should be 
included as well as changes in refining processes and capa- 
cities which have occurred since 1955* A model for JP5 
will be complicated by the fact that in time of war, other 
military fuels would be competing for a portion of the 
total product-mix and that industry and civilian require- 
ments would also be altered. Nevertheless it is considered 
that such a linear programming model can be constructed and 
manipulated to provide approximate answers to some of the 
questions which are now unresolved. The ability of the 
model to adjust for reduced capacities also permits con- 
sideration of alternative actions in the event of destruc- 
tion cf a portion of the nations refining capacity. 

On the basis of the author's review of information re- 
lating to the analysis of industry wide production and capa- 
city problems, it is considered that linear programming 
offers a generally satisfactory method of solution. Other 
attempts to predict the capabilities of an economy have been 
made using gross national product which usually lead to 
overly optimistic results when applied to short range prob- 
lems because they overstate the degree of substitutability 
of products. On the other hand, using a single area of pro- 
duction as a basis for prediction usually understates sub- 
stitutability and consequently under-estimates capabilities. 
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The Leontief interindustry flow model was an improvement 
and efforts to correct the errors introduced by rigid assump- 
tions in even this model have led to the use of linear pro- 
gramming as an extension of Leontief' s original methods. 
Linear programming allows both for the production of a 
specific item by any of several items by a single process. 

The development of an industry process analysis model 
is not a simple nor inexpensive task. Success requires 
close cooperation between the model builder, the sources of 
data and the numerical analyst. Perhaps the most critical 
area in constructing the model is in the collection of 
truly significant data and an understanding by the model - 
builder of the use of such data. While difficult, it is 
feasible to construct a model which will provide reasonably 
good results. The use and interpretation of the results 
then becomes the function of top management. 

Additional applications of linear programming in re- 
fining operation may be found in references 29 and j50- 



CHAPTER II 



GAME THEORY 
BACKGROUND 

Competitive situations are distinguished by the fact 
that two or more individuals are making decisions in situa- 
tions involving conflicting interests and in which the out- 
come is controlled by the decisions of all concerned. Such 
situations are found in many social, political, economic 
and military problems. Most of these situations involve 
elements of chance but in each case it may be assumed that 
each opponent will act in a rational manner and will attempt 
to resolve the conflict in his favor. 

Von Neumann developed an approach to competitive prob- 
lems utilizing the minimax principle which has as its basic 
idea the minimization of maximum loss. Von Neumann's ideas 
have led to the development of a branch of mathematics 
known as game theory in which the aim is to find the optimum 
strategy for a competition. 

The foundations of game theory were laid in the book, 
" Theory of Games and Economic Behavior 11 by Von Neumann and 
Morgenstem^ 1 , first published in 19^4. Newman^, in coment- 
ing on the social application of mathematics, wrote, - 
"the theory of games can fairly be said to have laid the 
foundation for the systematic and penetrating mathematical 
treatment of a vast range of problems in social science. " 
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However, a more recent survey of game theory by Luce and 
Raiffa"^ indicates that there have been few applications in 
the social sciences. The majority of applications have 
been in military operations and in competitive industrial 
problems such as bidding strategies. 

During its short history, game theory has had a num- 
ber of contributors. Among the more prominent are George 
W. Brown who developed a method of finding solutions by 
ficticious play, and J. D. Williams who defined a simple 
method of interpereting equations called "the method of 
oddments". H. Kuhn developed the "kernel method" for solu- 
tion of finite games and George Dantzig proved game problems 
could be converted to linear programming problems and vice 
versa. 

New mathematical tools are being created to develop 
and extend the applications of game theory in the social 
sciences, operations research, statistics and other fields. 

EVALUATION OP BIBLIOGRAPHY 

For those seeking references and articles in specific 
areas of game theory, the excellent bibliography by 
Batchelor-^ is recommended as well as the Index of Publica - 
tions of the Rand Corporation^ . The Rand Corporation has 
sponsored the work of a number of men in this area. Both 
sources are well indexed. 

Among texts dealing with the subject. Introduction to 

^36 

Operations Research by Churchman, Ackoff and Arnoff^ is 
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good and includes the most complete treatment of competitive 
bidding models encountered in this study. The book. Mathe - 
matical Methods of Operations Research by Saaty^^, is also 
recommended. 

A very different approach to game theory is found in 
Executive Decisions and Operations Research by Miller and 
Starr^. The treatment is rooted in decision theory and a 
general, non-mathematical, presentation is given. One of 
the better known works in game theory is the Compleat Stra - 
tegist by Williams^ which is suggested to those with a 
limited background in mathematics. 

A variety of articles on applications and extensions 

of game theory can be found in the pages of the J ournal of 
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the Operations Research Society of America and in Manage - 
4l 

ment Science . 

For those interested in applications of game theory in 

42 

the area of war games, the paper by Dresher , covers some 
interesting aspects of target selection and prediction as 
well as optimal military strategy and tactics. 



MATHEMATICAL METHODS 

In the theory of games, we are dealing with systems 
in which two or more decision makers are in a competitive 
situation and in which the outcome is controlled by the 
decisions of all players. While not all competitive situa- 
tions can be analyzed by means of game theory, there are 
many military, social, and economic situations which are 
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applicable. In line with the subject of this thesis, our 
purpose in examining game theory is to set forth methods 
for optimizing the strategy of a player. 

A competitive game has the following four properties: 

a. The number of players is finite. 

b. Each of the N players has a finite group of pos- 
sible courses of action which need not be the 
same for each player. 

c. A play in the game occurs when each competitor 
selects one of his alternatives. It is assumed 
that selections are made simultaneously and with- 
out knowledge of the choices of other players. 

d. The outcome of a play determines a set of pay- 
offs (positive, negative, or zero), one for each 
player. 

In a game where there are N players, and player i has 
n^ possible courses of action, there are . . . n^ pos- 



sible outcomes for each play. An outcome 0 results in a 
payoff R(i,0) to player i. If for every possible outcome 
0, we have 



then the game is called a zero-sum game. That is, if the 
sum of all payoffs to all players is zero, positive to some 
and negative to others, it is a zero-sum game. Poker is 
such a game. 



i=n 




( 2 - 1 ) 



i=l 
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The strategy of a player is the decision rule he uses 
for selecting a course of action. The. decision rule should 
not require definite information about an opponents choice. 

A pure strategy is a decision by a player in advance of all 
plays, to select a single course of action which he then 
follows for all plays regardless of the outcome. A mixed 
strategy is a decision, in advance of all plays, to select 
an alternative for each play in accordance with some parti- 
cular probability distribution. A mixed strategy is more 
advantageous than a pure strategy after the pattern of play 
becomes evident in that players are unable to predict the 
next play of their opponents. 

A pure strategy is frequently identified by a single 

number such as 1 or 2 representing the course of action 

selected from all available alternatives. A mixed strategy, 

for a player with m possible courses of action, is denoted 

by the set, X, of m non-negative numbers whose sum is unity, 

representing the probabilities with which each alternative 

is chosen. If x^ is the probability of choosing course i, 

then the set X is (x, , x 0 x ) 

12 m 

where x i 0 ( i = 1, 2, . . . . , m) 

i=m (2-2) 




i=l 



It will be recognized that a pure strategy is a special 
case of a mixed strategy where all of the x^ are zero except 
one which has a value of 1. A player may select any one of 



44 



m pure strategies but he has an infinite number of possible 
mixed strategies. 

Our purpose here is to select the optimal strategy 
for a player and for this we must have a criterion of opti- 
mality. The criterion normally used is known as the maxi- 
mum criterion. This can be explained in that every player 
must be regarded as desiring to maximize his returns and 
minimize his losses. These desires are in conflict with 
those of other players since some player must lose if 
another is to win. The solution then consists of specify- 
ing the strategies for each player which will maximize his 
return as a winner and minimize his loss if he is a loser. 

Two Person Zero-Sum Games 

A game with only two players A and B is most conven- 
iently described by means of a pair of matricies. Row 
designations for each matrix are the courses of action avail 
able to A; column designations are the courses available to 
B; cell entries are the corresponding payments to A for one 
matrix and to B for the other matrix . In a zero-sum two 
person game, the cell entry in B's payoff matrix is the nega 
tive of the corresponding entry in A's payoff matrix, and 
for this reason problems are usually worked out using only 
one matrix. 

All two person zero-sum games can be solved. The simp- 
lest type of game is one where the stable optimal strategies 
are pure strategies. This is the case only if the payoff 
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matrix contains a saddle point. The saddle point is the 
element of a matrix that is both the lowest element in its 
row and the highest element in its column. In other words, 
a saddle point is at once the largest of the row minima and 
the smallest of the column maxima. 

In such instances, the solution is for player A to 
follow a pure strategy corresponding to the row through the 
saddle point and for B to use his pure strategy correspond- 
ing to the column through the saddle point. The value of 
the game to a player is his expected gain in one play of the 
game when both players use their optimal strategies. 

For purposes of illustration, a payoff matrix is given 
below in which the row minima have been circled and the 
column maxima have been enclosed in, boxes. 



I 

II 

A III 
IV 




B 

III IV 

1 !"8 

6 

3 

2 2 





nn 

© 



It can be seen that a saddle point of 4 exists in this 
case. The optimal strategies of players A and B are thus II 
and III respectively. The value of the game to A is 4 and 
to B it is -4. When B learns that A is playing strategy II, 
he will follow strategy III because it minimizes his loss. 
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Mixed Strategies 

When there is no saddle point, the best strategies are 
mixed strategies and the solution to the game consists of 
evaluating the probabilities with which each pure strategy 
should be used. If a mixed strategy is utilized a minimum 
gain cannot be assured on any specific play, however, we 
attempt to insure that the minimum possible value of the 
expected gain will be as large as possible. 

Let us assume a competive two person zero-sum game 
with a payoff matrix as shown below. 



I 



I 



-1 



II 1 
III 3 



B 

II III 
2 1 

-2 2 

4 -3 



It will be noted that there is no saddle point. Search- 
ing for pure strategies, A finds I to be most advantageous 
and B selects III, assuring themselves of gains of at least 
-1 and -2 respectively. (it will be recalled that the payoff 
to B is the negative of that in A's matrix.) If the players 
start with their pure strategies, A soon discovers that B 
always chooses III and realizes that he. A, can increase his 
gain from 1 to 2 by playing II. B then shifts to II to in- 
crease his gain to 2 on each play. The strategy for each 
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thus becomes one of continuous change. It is also true 
that, if one player ever realizes that the other is using a 
fixed pattern or sequence of strategies, he will be able to 
anticipate the next play of his opponent. A good way to 
avoid a fixed pattern of play is to make the strategy for 
each play a random variable. 

Let X = (x^, x 2 , x-j) and Y = (y^, y^, y^) indicate the 
optimum mixed strategies for A and B respectively. The prob- 
ability that A will play i is x^ and the probability that B 
will play j is y^. The gain to A at each play now becomes 
a random variable a and the expected payoff of a play to A 
is: 

E(a;X,Y) = ^a (2-3) 

ij 

where a^ is the gain to A when A plays i and B plays j. 

Without regard to Y, A seeks to choose X so that his 
expected gain exceeds some quantity v^. The expected pay- 
off for B similarly becomes: 

E(-a;XY) (2-4) 

ij 

B plays so that his expected gain exceeds a quantity 
v 2 » Designating optimum strategies as X Q and Y Q for A and B 
respectively, we must then satisfy the following conditions: 



E(a;X Q ,Y) ;> v x for all Y 


(2-5) 


E(-a;X,Y Q ) £ v 2 for all X 


(2-6) 
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One of the theorems developed by Von Neumann states 
that In any zero-sum two person game, the quantities v^ 
and Vg above, are the negatives of each other. Substitut- 
ing v for v^ and -v for v^, the following relations are 
obtained. Note that in equation (2-8), reversing the in- 
equality sign causes the negative signs to be dropped. 



E(a;X Q ,Y) £ v for all Y (2-7) 

EfajXjY ) < v for all X (2-8) 

Therefore: 

E(a;X 0 ,Y 0 ) = v (2-9) 



The interpretation of this is that if both players 
use their optimal minimax strategies, their actual expected 
gains coincide with their minimax expected gains. 

Applying the foregoing to the numerical example on 
page and substituting the values of a., in (2-7 ) , we 

-L J 

may obtain an expression for the expected gain to A for 
this game. 

Y 1 (-x 1 + x 2 + 5x^) + y 2 (2x 1 - 2x 2 + 4x^) 

( 2 - 10 ) 

+ + 2x 2 - 3x^) ^ V 

Since the foregoing expression is true for all y ., it must 

<3 

be true when any two of the y^ are zero and the third is 
one. This observation allows us to create from the rela- 
tion in (2-10), three inequalities which are: 
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-x ± + x 2 + ^ v 

2x 1 - 2 x 2 + 4x^ ^ v (2-11) 

X 1 + 2x 2 " 5x 3 ^ V 



The first of these was obtained by setting y^ = y 2 = 0 
y^ =1. The second by setting y^ = y^ =0; y 2 = 1 and the 
third by setting y^ =1; y 2 = y J = 0. 

By a process similar to that used to obtain (2-10), 
we may write an expression for the expected gain to B and 
use this to write three more inequalities similar to those 
in (2-11) above. 

-y 1 + 2y 2 + y^ ^ v 

y x - 2J 2 + 2y^ < v (2-12) 

3y 1 + 4y 2 - < v 

It is also true that: 



x i ^ 0 


(i = 1, 


2, 3) 


<< 

-J. 

I(V 

0 


(i = 1, 


2, 3) 



Since the sum of the probabilities x^ is equal to 1 
from (2-2), we can also write: 



X 1 + x 2 + x 3 = 1 
y l + y 2 + y 3 = 1 



(2-14) 



The equations and inequations included in (2-11) through 
(2-14), describe the seven unknowns x^, y^ and v. It is 
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apparent that not all the and y^ can be zero, but it is 
possible for equality signs to hold in all the relation- 
ships of (2-11) and (2-12). This is first assumed to be 
the case and we can undertake a solution of the following 
set of eight equations with seven unknowns: 



-x 1 + + 3x^ = v 

2x 1 - 2x 2 + 4x^ = v 

X 1 + 2x 2 “ ^ x 5 = v 
~Y 1 + 2y 2 + y^ = v 

y l “ 2y 2 + 2y 3 = v 
3y 1 + 4y 2 - 3y 5 = v 



X 1 


+ 


X 2 


+ 


x 3 


= 1 


1 — 1 


+ 


y 2 


+ 


y 3 


= 1 



(2-15) 



If positive solutions are found for the x i and y^, 
then all relations are satisfied. If not, one or more of 
the equality signs in (2-15) must be replaced by inequali- 
ties until positive solutions are found. 

In this instance, the equations hold and positive solu 
tions are: 



x. 



L 1 2 
17/46 20/46 



x 5 

9/46 



y i y 2 y^ v 

14/46 12/46 20/46 30/46 



Note that A now has an expected gain per play of 30/46 
where with his best pure strategy he could only be sure of 
not losing more than 1. B now has an expected loss of 
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-50/46, where with a pure strategy he could not be sure of 
losing less than 2. 

Solution of Rectangular Games 

Finite two-person zero-sum games are referred to as 
rectangular games. Three theorems due to Von Neumann are 
basic to the solution of rectangular games and will be re- 
viewed here in brief. 

Theorem 1. Every rectangular game has a value and a 
player of a rectangular game always has an optimal strategy. 

Theorem 2. v*, X*, and Y* are respectively the value 

of a rectangular game and optimal strategies for players 
and Pg, if for pure strategies the following relations 
are true: 

E(Xp,Y*) < v* for every pure strategy Xp of P^ (2-16) 

E(X*,Y ) ^ v* for every pure strategy Y of P ? (2-17) 

Theorem 1 indicates that each rectangular game has a 
solution with a value equal to some finite number. The 
theorem does not say that each game has a unique solution 
and there may be more than one optimal strategy. 

Theorem 2 serves to check proposed solutions and if 
the relations do not hold for the proposed solution, it is 
not valid. 

With this background, let us examine briefly a few of 
the methods of solution which have been developed for rec- 
tangular games. 
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2X2 Games 

The general payoff matrix for a 2 X 2 game is shown 

here: 



P 



2 



Alternative 




B 2 


i a 2 

I 


a ll 


a 12 


j A 2 


a 21 


a 22 



If the pay-off matrix has a saddle point , the solution 
is direct and need not be considered further. However for 
purposes of discussion let us assume that the game has no 
saddle point. For this situation it may be shown that the 
specific relationships among the elements a^ . are: 



a ll 

a 22 



> a 12 


a ll >a 2l 


> a 12 


a 22 > a 2l 



(2-18) 



If these relations were not true, a saddle point would 
exist. 

Let the general solution of this game be 

V X o ■ 1_x o ) ' Y o = 1_y o ) 

where 0 £ x Q £ 1 

and 0 £ y Q £ 1 



( 2 - 19 ) 
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By theorem 2: 

a ll y o + a 12 (1 ' y o* S v o for A 1 

a 21 y o + a 22 (1 - y o> S v o for A 2 

(2-20) 

a ll x o + a 21 i v o f0r B 1 

a 12 x o + a 22 (1 ' x o> i v o for B 2 

It* 

It has been demonstrated ^ that these four inequali- 
ties have a solution only when the equalities hold, the solu- 
tion being: 

Vq = [ ]_®'22 ” a 12 a 21 ^ ^ a ll a 22^ ” ( a i2 a 21^“ 

( 2 - 21 ) 

X o = ^ a 22 “ a 2l)/^ a ll + a 22^ " ^ a 12 + a 21^ 

( 2 - 22 ) 

y o = ^ a 22 " a 12^ // ^ a ll + a 22^ " ^ a 12 + a 21^ 

(2-25) 

This means that players and P 2 should use the following 
optimal mixed strategies: 

x _ r a 22 “ a 21 

0 (a ll + a 22 J - la 12 + a 21> ’ 

( a ll + a 22^ ' (a 12 + a 21^ f ° r Pl 



(2-24) 



5 ^ 



Y _ r a 22 " a 12 

O 4a n + a 22 ) - (a 12 + a 21 ) 



( 2 - 25 ) 



(a n + ^J !la lg Va^T 1 fOT P 2 



A method for interpreting these two equations was pre- 

iiii 

sented hy J. D. Williams , which he terms the method of 
oddments. Williams calls the difference between two num- 
bers in a row or column an oddment. He terms the ratio 
of the two relative frequencies of any player, odds. The 
optimal mixed strategy for thus is to play A 1 and Ag 
with the odds of oddment a 22 - a 21 to oddment a^ - a ]_2* 

The optimal mixed strategy for player P 2 is to play and 
Bg with odds of oddment a 22 - a^ 2 to oddment a^ - a^. 

In more familiar numerical terms, if the payoff matrix is 
such that: 

a 22 " a 21 = 5 

a ll " a 12 = 5 

Then the optimal odds for player P^ are 5 to 5 and he 
should thus play alternate A 1 3/8 of the time and A g 5/8 
of the time. 



M X N Games 

Several methods have been developed for solution of 
games of this type. All are relatively long and no attempt 
will be made to cover them in detail here. A few will be 
briefly described. 
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In selecting an approach to a game and before attempt- 
ing a solution, two steps should be taken. As stated earlier, 
the first thing to do is to examine the problem for a saddle 
point. Failing this, certain alternatives can be eliminated 
by checking for what is termed dominance. When the payoffs 
for alternative are greater or equal to each of the cor- 
responding payoffs for Ag, A^ is then said to dominate and 
Ag may be disregarded. The concept of dominance may be 

extended to more than two alternatives at a time. McKinsey 
lifi 

and Von Neumann give detailed examples of this approach. 

Ficticious Play Solution 

b7 

George W. Brown's 1 method of solution by ficticious 

play is based on a hypothetical series of consecutive plays- 

of the game and uses an infinite sequence of steps to reach 

a solution. In the first play, P-^ selects A^ and Pg plays 

B-^. Thereafter each player chooses the optimum pure strategy 

against the mixture represented by all of his opponents 

48 

past plays. Saaty presents a good example of this method. 
There are other techniques of solution including the geomet- 
ric, algebraic, matrix and Kernel methods. For examples and 

for further information on these- methods, the reader is re- 

49 SO 

ferred to Saaty ^ or McKinsey-^ . 

Linear Programming Solution 

It has been shown by Dantzig-^ that every game problem 

can be transformed into a linear programming problem and 

S2 

vice versa. Saaty-' presents some examples of how this is 
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accomplished and it is sufficient for us to note that 
Dantzig has proven the validity of this transformation and 
it becomes a valuable tool in both the solution of game 
and linear programming problems. 

N Person Zero-Sum Games 

It should be noted that we have been dealing with rela- 
tively simple games involving two players. No generally 
acceptable definition of a solution of an n person zero- 
sum game yet exists. One reason for the lack of an accept- 
able definition is that coalitions can form when more than 
two players participate. Von Neumann and Morgenstern^ 
have suggested an approach based on possible formation of 
coalitions involving dominance. However their definition 
of a solution does not show how to find optimal strategies 
nor indicate what the final outcome will be. For these 
reasons Von Neumanns definition of n-person games has not 
had wide acceptance. 

Non-Zero Sum Games 

A discussion of the methods of game theory would not 
be complete without some mention of non-zero sum games. 

Here too, there is no generally acceptable definition of a 

54 

solution. Von Neumann suggests an approach introducing a 
ficticious player assigning him such payoffs that the game 
is reduced to a zero-sum game. The solution is then 
attempted by methods applicable for zero-sum games. However 
there are infinitely many solutions and no acceptable defini 



tion has been found. 
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APPLICATION OP GAME THEORY 

The most extensive applications of game theory have 
been in the area of military warfare. It has also been 
applied to economic problems involving commercial compe- 
tition and in both cases the complexities of the problem 
must be simplified by various assumptions to make it amen- 
able to analysis. This creates some doubt as to the valid- 
ity of the results, especially where the game gives nearly 
equal opportunity to the players. In the case of war games, 
frequently one side has a considerable advantage and the use 
of game theory analysis will usually reveal this advantage 
which might otherwise remain obscured. 

In the area of military petroleum problems, the petro- 
leum logistics involved would usually be a portion of the 
battle game which will not be treated here. There is, how- 
ever, an application of game theory in the area of compe- 
titive bidding which has a direct bearing on military petro- 
leum procurement which will be presented here as an example. 
While this competitive situation differs from the strict 
sense of a competitive game, it is possible to use similar 
means of analysis. 

Competitive Bidding 

As a simple case, let us assume a refiner engaged in 
bidding on military petroleum contracts, makes it a practise 
to record, in addition to his own bid, his cost estimate K 



and the lowest bid B by his opponent. The refiner finds 
that the random variable x, where 
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follows a normal distribution f(x) with a mean of 0. 15 and 
standard deviation of 0.05. 




Figure 5 

Graphical Solution of Equation ( 2 - 32 ) 



If the refiner is interested in maximizing his expected 
profit on each contract, and if it is assumed that his com- 
petitors bidding habits will not change, what should the 
refiners bidding strategy be? 
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Let A be the refiners 1 bid and B be the lowest compet- 
ing bid. The refiners' profit P is thus: 

P = A - K if A < B 

(2-26) 

P = 0 if A > B 

Let E(P;A) signify the expected profit with bid A. 

Then: 

E(P; A) = (A - K)P if A < B (2-27) 

Replacing the decision variable A by an equivalent variable 



Z, where 




CS3 

II 

> 
^ 1 


(2-28) 


Substituting from (2-28) 




E(P; A) = KZP if Z < x 


(2-29) 


r 00 

E(P; A) = KZ J f (x) dx 


(2-50) 


Z 





Where x is the random variable and f(x) is its function. 
The quantity . 

00 

Z J f (x ) dx (2-51) 

Z ‘ 

is maximized by the solution Z of the equation 
n 00 

J f (x ) dx - Z f(Z) = 0 (2-52) 

Z 

Graphically, E(P; Z) is a maximum when the two cross 
hatched areas in Figure 5 are equal. 
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Numerically, Z may be found by using tables of ordinates 
and areas of normal distribution. In this instance, 

Z = 0.117 is a solution to the foregoing equation. Thus 
the refiners optimum bidding strategy is to add 11.7$ to 
his own cost K, in preparing his bid. 

A critical assumption is that the competitors will not 
change their bidding habits. In the case where only two 
bidders are in competition, such an assumption will not 
hold, however, it is nore likely that where several competi- 
tors exist, it may be more reasonable to assume that the 
bidding pattern will not change. 

In the case of petroleum refiners, competing companies 
have a fair amount of information on their opponents capa- 
bilities and markets plus a reasonably good estimate of 
their costs. This information enables them to bid more 
intelligently, no doubt, than the refiner in the example. 

It is known that petroleum refiners in at least a fev; cases, 
have prepared programs for computer calculation of bids. 

It is not surprising that these computer programs are not 
available for analysis. 



CHAPTER III 



INVENTORY THEORY 
BACKGROUND 

Probably more work has been directed toward inventory 
control than toward any other problem area in business, 
industry, and government. In spite of this, or perhaps 
because of it, a problem in terminology exists in the area. 
Moreover, there are widely divergent views on what research 
in this area should include. Some references treat inven- 
tory systems as record-keeping problems involving quantities 
of stock while others take a broader view and consider the 
overall problem of stocks and inventory investment. Still 
others are concerned with what items to stock, quantities, 
and utilization of facilities. The diversity stems from 
the many facets of what has been called the inventory sys- 
tem and the large variety of systems found in practise. 

In this thesis the basic characteristics of inventory 
systems, input, inventory levels, and output as well as the 
costs associated with these factors, will be discussed and 
and applied in a few models. The aim in the handling of 
each inventory problem is one of optimizing costs within 
the limits of effectiveness required by management policies. 

Historically, one of the earliest efforts at analysis 
of inventory systems was made by F. W. Harris in 1915. 

Harris recognized many of the elements of inventory problems 
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and developed a formula for the determination of optimal 
lot size. ^ 

q o = ^ 2rc 5 /c 1 (5-1) 

In this formula q is the economic lot size, r is the rate 

o 

of requirements per unit time, c^ is the carrying cost of 
one unit in inventory per unit time, and c^ is the cost of 
replacing inventory. Application of this formula has been 
extensive; however, Harris assumed that the rate of require- 
ments were known, that no shortages occurred, and that lead 
time could be disregarded. These limitations made it use- 
ful only in special cases and improvements were sought. 

Cooper in 1926 and Pry in 1928 analyzed systems with 
varying production and demand. A significant advance was 
made by Pry when he showed that the theory of probability 
can be applied to some inventory systems when requirements 
are not known with certainty. 

Another important step was taken by Raymond who made 
the first effort to analyze and develop a theory of inven- 
tory systems in his book Quantity and Economy in Manufac - 
ture published in 1951 and no longer in print. 

Dvoretsky, Kiefer and Wolfowitz^ published a paper, 
"The Inventory Problem", in 1952 in which they studied 
systems of a general nature utilizing extensive mathemati- 
cal and statistical methods and endeavored to provide solu- 
tions for general cases. While they were successful in 
obtaining some unique solutions, there is some danger in 
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applying general solutions to the analysis of specific 
inventory systems and their work is chiefly of theoreti- 
cal interest. 

s8 

The approach taken by Arrow, Harris, and Marschak^ , as 
well as Whitin'^, was to analyze specific hypothetical sys- 
tems encountered in practise and apply the explicit rules 
developed to actual cases. 

A different approach by Naddor and others, places 
emphasis on the development of methods of analysis. Their 
intent is to develop models of elementary systems and analyze 
more complicated systems as extensions of the basic models. 

In reviewing the numerous sources of information on inventory 
theory, the text by Churchman, Ackoff, and Arnoff° is con- 
sidered the best for those seeking basic knowledge on in- 
ventory and purchasing. Another recommended source is Pro- 

62 

duction Planning and Inventory Control by Magee 

A review of the references and bibliography at the end 
of this thesis will suggest additional sources of informa- 
tion. The periodical. Management Science , is an excellent 
source of information on inventory theory and other related 
subjects. 

ANALYSIS OP INVENTORY SYSTEMS 
An inventory system is usually considered to be one in 
which the following types of costs are significant, and in 
which any two or all three are subject to control: 




















( 
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1. the cost of carrying inventories, 

2. the cost of incurring shortages, 

J>. the cost of replenishing inventories. 

The first cost relates to inventory investment, stor- 
age and handling costs, and losses due to spoilage, evapora- 
tion, etc. The second, shortage cost, concerns the cost due 
to lost sales, overtime, special handling, etc. The third 
cost is that of ordering, machine set up, etc. The general 
terms for the three costs are surplus cost, shortage cost 
and set-up cost. As will be noted in some applications in- 
volving use of inventory theory, all other costs become 
significant in actual cases and problems do not limit them- 
selves to the area of inventories alone. Actual problems 
frequently require use of other optimization techniques in 
order to avoid what we have termed sub-optimization. 

Inventory Problems 

An inventory problem is one of making optimal decisions 
with respect to an inventory system in order to minimize the 
total cost of the system. The measure of effectiveness in 
inventory problems is one of costs. It is assumed that all 
shortages, surpluses and set-ups can be expressed in common 
units of the measure of effectiveness even though the units 
may not be dollars. 

Control of Variables 

While it has been indicated that the surplus, shortage 
and set-up costs are the factors which we desire to control. 
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the optimal decisions designed to minimize total costs are 
normally made in terms of time and quantity. More specific- 
ally, the decision is frequently one of when to replenish 
the inventory and by how much. The variables thus subject 
to control are usually time and quantity. 

Elements of Inventory Systems 

The three elements of inventory systems are output, 
input, and costs. Output refers to what is removed or 
issued from inventory. Input includes what is received or 
taken into inventory and costs refers to the three costs 
mentioned earlier (surplus, shortage and set-up costs). 

Output : In several respects, the output of a system is the 

most important element. Inventories are maintained to fill 
demands. The demand or output cannot usually be controlled 
but the pattern of output in terms of time and quantity 
has a significant bearing on the operation of the system. 

The output may be variable or constant. When accurate ad- 
vance information about the size of the output is available, 
it is said to be known with certainty. When not known it 
may be subject to a known probability distribution. The 
output rate or demand rate is the output size per unit 
time and it may be constant or vary in a manner similar to 
the output size. 

Output patterns vary widely but may generally be de- 
scribed mathematically if sufficient data are available. If 
the output size is x during an interval of time t, the output 
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pattern may follow one of the patterns in Figure 6. All x 
units may be drawn at the beginning or at the end of the 
period, uniformly during the period, or in accordance with 
some other distribution. 




Figure 6 
Output Patterns 



The output patterns shown can be described by: 

n 

q(t) = s - x nT r/t ( 3 - 2 ) 

where Q(T) is the total inventory at time T 
S is the inventory when T = 0 
x is the output during the interval t 
n is the output pattern index 
t is the scheduling period 
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Input : The input element can generally be controlled. In- 

put involves not only the quantity scheduled to be taken 
into inventory, but the time of ordering or scheduling as 
well as the time of receipt into stock. The scheduling 
period is the term used for the period between decisions 

r 

regarding inputs. This cannot always be controlled. Lead 
time is the period between scheduling or ordering and re- 
ceipt and this too is frequently not subject to control. 

The input size is defined as the quantity ordered or 
scheduled for input. Numerous inventory problems are con- 
cerned with the optimization of the input size in which 
case it is frequently called the optimal or economic lot 
size. 

Input patterns are related to the input period which 
is the time interval t' during which the input size q is 
being added to the inventory. The average input rate p 
is the ratio of input size to input period or units of in- 
put per unit time. 



Input patterns vary widely as illustrated in Figure 7 
and each type is important in some particular problem. 

All of the patterns shown have an input size q Ijeing 
added at average input rate p. This group of patterns can 
be described by: 



p = q/t' 



(3-3) 




(3-4) 
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where: Q(T) Inventory at time T 

q Input size 
t ' Input period 

p 'Average input rate (p = q/t 1 ) 
n Input pattern index 




Figure 7 
Input Patterns 



Costs : The cost element of inventory systems is considered 

to have been covered sufficiently in the section headed Analy- 
sis of Inventory Systems on page 63 - If any further clarifi- 
cation of the types of costs included under the titles of 
surplus costs, shortage costs, and set-up costs is desired, 
the following references are recommended: 64, 65 . 
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APPLICATIONS OP INVENTORY THEORY 

In order to illustrate the application of inventory 
models, a few types will be presented here, largely in the 
form of examples. Among basic inventory systems, there 
are three types all of which involve the balancing of two 
of the three types of costs; surplus, shortage and set-up. 
Balancing Surplus and Shortage Cost s 

The following conditions are known: 
c^ (set-up cost) = constant 
t (scheduling period) = constant 

IT 

r (output ratio in units per unit time) = constant 
c^ (surplus cost of one unit per unit time) = constant 
Cg (shortage cost of one unit per unit time) = constant 

The inventory level S at the beginning of each schedul- 
ing period is the only variable subject to control in this 
model. The system may be illustrated graphically as shown 
in Figure 8. 




System with Surplus and Shortage 
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The sawtooth line in the figure represents the inventory 
level at various times. The initial inventory level S is 
drawn down at a constant rate ir. At time the inventory 

reaches zero and a shortage s occurs. The demands which 
represent shortages are accumulated and stock is issued to 
cover these demands after the inventory is replenished. In- 
put q is sufficient to cover the shortage and raise the in-, 
ventory level again to S. Note that the vertical input 
line indicates practically instantaneous receipt of stock 
at each scheduling period t . 

Since shortages are carried over as indicated, the 
input q becomes: • 

q = rt p (3-4) 

where r is the rate of output and t p is the scheduling 
period or time betxveen inputs. 

The inventory level, s, which is actually the amount 
of shortage at the end of the scheduling period, is depen- 
dent upon the variable S as well as on the other conditions 
given. Although s is a negative inventory, it is treated 
as a positive quantity here: 

s = rt p - S (3-5) 

The solution of this system requires the determination 
of S Q , the optimal inventory level, and C , the optimal 
(minimum) total cost per unit time. 

The costs of the system are defined by the equations: 

(3-6) 



Surplus cost per unit time 
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Shortage cost per unit time = IgC^ (3-7) 

Setup cost per unit time = Nc^ = c^/t (3-8) 

where 1^ and I 2 respectively, are the average surplus and 
shortage during a unit of time and N is the number of set- 
ups per unit time. 

Using equations (3-8), (3-7)# and (3-8) as well as 
the relationships from Figure 8, one may derive the follow- 
ing equations which describe the surplus costs and shortage 
costs of the system: 

C 1 = (l/2)S(t 1 /t p )c 1 = S 2 C] /(2rt p ) (3-9) 

C 2 = (l/2)s(t 2 /t p )c 2 = s 2 C2 /(2rt p ) (3-10) 

The derivation of (3-9) can be explained briefly by 
pointing out that the term (s/2) represents the average 
surplus inventory and that the term (t^/t ) represents 
that proportion of the scheduling period during which a 
surplus exists. The product of these terms multiplied by 
c^, the surplus cost of one unit per unit time, thus equals 
the total cost of surplus C^. The alternate form of the 
equation is obtained by substituting the relation t^ = S/r 
which is the time required to deplete the initial inventory 
S at. output rate r. Equation (3-10) is derived similarly. 

Adding (3-9) and (3-10) one may obtain the total cost 
per unit time C which is subject to control since C^, the 
setup cost, is constant. 

C = C 1+ C 2 = [S 2 c 1 + (rt p -S) 2 c 2 J/(2rt p ) 

= [l/(2rt p )]( C;L +c 2 )S 2 - c 2 S + Cgrtp/2 (3-11 ) 
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This equation, usually called the total cost equation, 
can be used to arrive at the optimal decision rules by dif- 
ferentiating, setting the derivative dC/dS equal to zero, 
and solving for S , the optimal inventory level. The opti- 
mal inventory level then becomes: 



This value of S can be substituted in equation (3-11) to 

obtain the minimum total cost, c . 

7 o 



Pump Seal Problem 

The foregoing model will be applied here in an example 
to illustrate the use of the equations. 

A fuel terminal uses 2400 pump seals annually due to 
rapid wear. When a seal is needed but not available, leak- 
age occurs costing $9 per month. 200 seals are ordered each 
month and delivered in one day. Each seal costs $30. In- 
ventory carrying cost is 15 $> per year. 

In terms of the designations given: 

(Scheduling period) t = 1/12 year 

Sr 

(Output rate) r = 2400 per year 

(Unit surplus cost) c^ = $80 x 0.15 = $12 per seal 



S o “ rt p c 2^ c 1 +c 2^ 



(5-12) 



c 0 = (l/2)(rt p c 1 c 2 )/(c 1 +c 2 ) 



(5-13) 
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Prom equation (3-11 ) 

c = 0.5 S 2 - 108 S + 10,800 

If no shortages were permitted, S would have to equal 200 
seals in which case the total cost would be $1200 per year. 
However, when shortages are allowed and surplus and shortage 
costs are balanced, the optimal inventory level, from equa- 
tion (3-12), becomes; S Q = 180 seals. The minimum total 
cost for the year is then $1080 per year. 

Balancing Surplus and Set-up Costs 

An inventory system in which the cost and consequences 
of shortage is so high in terms of lost operating time, etc. 
that shortages are unacceptable, - may still be optimized by 
balancing the surplus and setup costs. Such a system is un- 
doubtedly appropriate for most packaged petroleum products 
for military use in that a shortage of a relatively inexpen- 
sive special purpose grease may prevent the operation of an 
aircraft or other vital equipment. In setting up this parti 
cular model, no reserve inventory level is included. This 
would not be the practise in an actual case but since a 
reserve level would simply constitute the addition of a 
constant cost, the treatment is no less applicable. An 
illustration of such a system is shown in Figure 9. 



Inventory 



7 ^ 




Figure 9 

System with Surplus and Setup Cost 



The conditions describing the system are: 

(No shortages allowed) = 0 

(Output rate in units 
per unit time) r = constant 

(Cost of carrying one 
unit per unit time) c^ = constant 

(Cost of each setup) c^ is known 

The variables in this model are the input q and the 
scheduling period t^. The relationship between the vari- 
ables is: 

q = rt p (3-14) 

It is assumed that the variable subject to control is 
q the input quantity. In this system, since no reserve 
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exists, the initial inventory S is equal to q and the short- 
age carryover, s, is zero since no shortage is allowed. 

The average inventory level during the scheduling 
period t is 1^ = q/2, thus by. equation (3-6) 

C 1 = qc x /2 (3-15) 

and by equations (3-8) and (3-1^) one may write: 

c 5 = CjAp = rc ^A (3-16) 

The total cost equation then becomes: 

C = C 1 + = qc 1 /2 + rc^/q (3-17) 

By differentiating and setting the derivative dC/dq equal 
to zero, we can solve for -the optimal value of q which gives 
the minimum total cost C. The optimal input value, q, is 
thus: 

q = */ 2rc^/c 1 (3-18) 

Equation (3-18) is frequently termed the "economic lot size" 
formula. 

Entering q in (3-1^) and (3-17) is then: 

t Q = ^ 2c^/(rc ;L ) (3-19) 

C o = J 2rc 1 c 5 (3-20) 

The two equations (3-18) and (3-19), provide the opti- 
mal lot size and optimal scheduling period as a basis for 
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the decision rules in operating the system in which surplus 
and setup costs are balanced while no shortages are allowed. 
Equation (3-20) describes the minimum cost of such a system. 

Pump Seal Problem (No shortage allowed) 

The model of the foregoing system, in which no short- 
ages were allowed, can be applied in the pump seal problem 
utilized in the previous application. In this instance we 
will assume that the scheduling period, t , is variable and 
not fixed; that no shortages are permitted; and that the 
terminal operator decides to make the seals with a setup 
cost of $350 for each production run. 

In this case it is necessary to decide how many seals 
should be made in each production run and how often this is 
to be done. 

Prom the previous problem, it is given that: 

Output rate r = 2400 seals per year 
Surplus cost c^ = $12 per seal per year 
Setup cost c^ = $350 per setup 

Substituting in equations (3-l8), (3-19), and (3-20): 

Optimal lot size = q = 374 parts per run 

Optimal scheduling 

period = t = I .87 months 

Optimal total cost = C Q = $4,488 

Application Involving Probabalistic Output 

An inventory system with a probability distribution of 
output is illustrated graphically in Figure 10. 
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System with Probability Distribution of Output 



Note that during each scheduling period, tp, the system 
may either have an output size x which is smaller or equal 
to the initial inventory level S, or it may be that x is 
larger than S. For each of these two cases the distribu- 
tion appears as shown in Figure 11. 




Figure 11 

Typical Cases with Probability of Output 



73 



< 

the pump seal problem, the following probability distribu- 
tion is assumed: 



Output size during 
Scheduling Period t 



Probability of 
Occurrence 



Probability Distribution of Output 

In this problem it will also be assumed that the 
scheduling period, t , is fixed at one month. It is also 
given that: 

Unit surplus cost c^ = $12 per unit per year 

Unit shortage cost c^ = $108 per unit per year 

Calculations are presented here in tabular form. The 
descriptions entered at the head of each column will assist 
in following the computation. 

The computations were initiated by assuming various in- 
ventory levels, S, and values of output, x, varying from 0 
to 500. One can then calculate the average inventory for 
cases where there is no shortage quite simply. For example, 
when S is 100 and x is 100, the inventory is drawn do wn to 
zero at the end of the scheduling period and the average 
inventory for the month is 100/2 = 50* When there is a 
shortage, as when S is 100 and x is 200, the inventory is 
drawn down to zero in the first half of the month. In this 



X 


0 


100 


200 


300 


400 


500 


P(x) 


.15 


. 20 


• 35 


.15 


.10 


.05 



Table 1 
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Inventory Level 


Output Size 


Probability of 
Output 


Average Surplus 
Case A 


Average Surplus 
Case B 


Average Shortage 
Case B 


Surplus Cost 


Shortage Cost 


Total Expected 
Cost 


Equation (3-21) 
Equation (5-22) 


Equation (3-23) 
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.20 


200 






2400 




480 








200 


• 35 


150 






1800 




650 


1560 




- 


500 


.15 




104 


4 


1248 


452 


252 








400 


.10 




78 


28 


938 


3024 


396 








500 


.05 




62 


62 


744 


6696 


372 


1020 


2580 


500 


0 


.15 500 






5600 




540 








100 


. 20 


.250 






5000 




600 








200 


.35 


200 






2400 




840 








500 


.15 


150 






1800 




270 


2250 






400 


. 10 




112 


12 


1544 


1296 


264 








500 


.05 




90 


40 


1080 


4520 


270 


534 


2784 

1 



Table 2 



Solution of System with Probability Distribution of Output 
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One may assume that an output of size x will occur dur- 
ing t with a probability P(x). When x £ S the average sur- 
P 

plus inventory is S-x/2 and the total surplus cost per unit 
time is (S-x/2)c^. The total expected cost C & for Case A 
(where all values of x £,S) is thus: 



x=S 



c a = C ]Z ( s_x / 2 ) 



(5-21) 



x=x n 



In Case B, when x > S, the total expected cost C, 
from Equation (5-11) becomes: 

x=x max • • 

C b = £ (IS\ + (x-S) 2 c 2 ]/(2x) | P(x)(5-22) 
x=Stl ^ ^ 



Note that the output size x = rt . In this derivation x is 

r 

substituted for (rt^) in Equation (3-11 )• 

The total cost equation then becomes: 



x=S 



x=x. 



max 



C = 



C o +0, 

a b 



= C 1 Z ( s_x /2)P(x) + c 1 ^ (l/2) S 2 [P(x) ] 



x=x 



x=S+l 



x=x max 

' + c 2 £ (l/2) (x-S) 2 [p(x)/x ] * (3-23) 

x=S+l 



Pump Seal Problem with Probability Distribution of Output 
In actual cases, the output of an inventory system is 
rarely constant but if sufficient data are available, a prob- 
ability distribution for output size can be determined. In 
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case the average inventory is 100/2 = 50 for only one half 
of the month so for purposes of this computation, the actual 
average during the scheduling period of one month is 25* 
Columns C 1 and C 2 are obtained by multiplying the average 
surplus or shortage by the unit surplus or unit shortage 
costs as applicable. The remainder of the calculation is 
generally straight-forward. 

In order to obtain the optimum inventory level S Q for 
this problem as well as the corresponding minimum cost, it 
is necessary to use interpolation and some trial and error 
calculations. In this instance, the approximate values of 
S Q can be seen to lie between 200 and 300 and comparatively 
close to 250. By additional computation it appears that S Q 
is approximately 235 seals at which inventory level the 
minimum total cost C Q is $2576 per year. 

Application of Inventory Theory in Military Tanker and Termi - 
nal Operations 

One of the important problems facing military and in- 
dustrial petroleum management, is that of optimizing the 
operation of systems encompassing the movement, storage and 
distribution of bulk petroleum products. The large amount 
of capital and effort expended in the operation of such 
systems, makes them an. important field for analysis. 

The military services are concerned with systems in- 
volving the transportation of products by tanker from re- 
fineries to terminals for storage and distribution. The 
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petroleum industry operates similar systems in that refined 
products are transported by tanker to regional terminals 
for storage and distribution. The industry operates another 
such system in the tanker movement of crude oil from 
terminals to refineries for processing. As operated by 
industry, these two systems are, of course, related in that 
the output of the crude oil system becomes the input of the 
refined product system. 

The similarity of the military and industrial refinery - 
tanker - terminal systems, suggests that optimization tech- 
niques applied by industry may prove valuable to the mili- 
tary services. '.It is the purpose and intent of this thesis 
to present such applications, and we will proceed on the 
assumption that recognition of the basic factors involved 
will assist toward a solution. 

Inventory - Queueing Systems 

Tanker-terminal systems are essentially inventory sys- 
tems, however, the problems associated with the arrival of 
tankers for service at a refinery or terminal involves cer- 
tain aspects of queueing theory. Since the subject of queue- 
ing has not been covered in this thesis, it is considered 
that a brief discussion of the elements of queueing situations 
and their solution by Monte Carlo methods should be included 
at this point as an aid to understanding what will follow. 
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Elements of Queueing Theory 

Queueing or waiting line theory concerns the attempt 
to describe and handle problems of organization and planning 
in the face of randomly fluctuating demands for a servive 
being performed. Situations frequently occur in which units 
arriving for service must wait before being serviced. If 
the conditions controlling arrivals, service times, and the 
order of servicing are known, then the nature of the waiting 
line can be analyzed mathematically. 

Two general classifications can be made, one a determi- 
nate type of operation and the other indeterminate. In the 
former all aspects of the system are known as functions of 
time and the state of the system can be predicted at any 
time. In the indeterminate case, aspects of the system are 
known only in terms of probabilities. 

The several characteristics of a queue, such as the num- 
ber in line at any instant or the waiting time experienced 
by a particular arrival, are random variables rather than 
being functionally dependent on time. Arrivals and service 
times are also considered to be random variables. Thus in 
determining the values of the characteristics of a queue, 
one is concerned with estimating only the average number of 
units in the line at any instant, the average time the ser- 
vice facility is idle during a day, the average service time 
and the average waiting time. 

In attempting to optimize a queueing system, the deci- 
sion maker uses his knowledge of the average characteristics 
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of the queue. He may consider: altering the number of ser- 
vice units, changing the average service time, etc. 

Poisson Arrivals 

For the most part, arrivals do not occur at regular 
time intervals but tend to be grouped or scattered in an 
uneven manner. Individual Poisson arrivals are completely 
independent one another. The assumption of Poisson be- 
havior of arrivals requires the presence of a constant X 

I 

which is independent of time, length of queue, or other 
random characteristics of the queue. The constant X repre- 
sents the rate at which units arrive for service. 

Monte Carlo Methods in Queueing 

Monte Carlo methods are very useful in queueing prob- 
lems which are difficult to analyze mathematically. The 
Monte Carlo method uses random sampling to play a game with 
a system in which an experiment is simulated. In the prob- 
lem at hand, the technique consists of simulating the arri- 
val of tankers by the use of random sampling. This method 
of simulated sampling has some advantages over actual samp- 
ling of queue characteristics, in that when done on a digi- 
tal computer, years of data can be developed in a few minutes 
and those factors which are subject to control can be manipu- 
lated. The effects of adding new piers or new tanks to a 
terminal system can thus be assessed on paper without actual 
installation or disruption of services. Saaty^ gives an 
example of a Monte Carlo solution in a queueing problem. 
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Limitations of Queueing Theory 

Queueing models should be regarded as approximations 
which serve primarily as guide lines. Arrivals in a sys- 
tem may vary from Poisson loading, which is frequently 
assumed, but a Poisson input can be regarded as the limit 
of disintegration for a scheduled process and thus pro- 
vides a backstop to keep the solution within bounds. An 
attempt should always be made to solve the queueing problem 
at its source, the random input. Any action which can be 
taken to schedule or control the input will be helpful. 

With this brief background in queueing theory and the 
Monte Carlo method, we can proceed to examine the problems 
of tanker- terminal systems. 

Tanker-Terminal Systems 

The similarity between military and industrial tanker- 
terminal systems becomes apparent when it is noted that each 
must maintain sufficient working or shipping tankage and 
pier facilities to handle the loading and discharge of tankers. 
The refiner maintains such facilities for the input of crude 
oil as well as additional tankage and pier facilities for 
the output of refined products. Moreover, the refiner must 
maintain additional tankage for seasonal storage of products 
such as heating oil. Tankage, comparable to the seasonal 
tankage operated by the refiner, is operated by the military 
in that extensive tankage is maintained at military terminals 
as reserve storage for emergency use. 
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The similarity of tanker and barge lifting of products 
from military terminals and refineries is less notable but 
still generally comparable for purposes of this application. 

In the servicing (loading and discharge) of tankers, where 
pier and pumping facilities are limited, delays may occur 
which are characteristic of queueing situations. The 
arrival and servicing rates of tankers becomes a signifi- 
cant part of the problem. In most cases, the queueing situ- 
ation is not believed acute at military terminals, but the 
military system is affected in that tankers arriving at 
refineries to lift military products become a part of the 
queueing problem at the refinery thus having some effect 
on arrivals at military terminals. In times of emergency 
this problem could become significant. 

Study on Shipping Tankage 

Although some attention has been given to the schedul- 
ing of tanker movements as transportation and linear pro- 
gramming problems, the published articles on tanker-terminal 

6 ? 

systems are extremely limited. An article by J. C. Dickson 
describes an operations research effort aimed at determining 
the optimal volume of shipping or working tankage at a re- 
finery to handle the loading of tankers. The system described 
includes the operation of additional seasonal tankage which 
is normally used for holding heating oil in storage to meet 
seasonal demands. The system treated _by Dickson is considered 
sufficiently similar to the military system to illustrate the 
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approach to the problem and the basic factors involved. 

The model developed by Dickson has been altered somewhat 
in order to fit the military system more closely and is 
discussed here. 

Description of the Problem 

A 

Since tankers do not arrive at fixed intervals and do 
not load or discharge a fixed quantity of petroleum pro- 
ducts, terminals and refineries must maintain shipping or 
working tanks to absorb the fluctuations experienced in 
input and output. The consequences of insufficient tankage 
may be delay of tankers, inability to load or accept a 
tanker, etc. In the case of the refiner, it may also mean 
the necessity for altering the entire production schedule 
of the refinery. 

The problem as seen by Dickson was one of minimizing 
costs. The costs considered to be pertinent to the problem 
were: 

1. Cost of delaying tankers 

2. Operating cost of working tanks plus carrying 
costs on these tanks 

3. Cost of reducing input 

4. Cost of pumping product to or from reserve storage 

5. Carrying costs on product held in inventory 

A schematic representation of the system altered to 
correspond to the military petroleum terminal operation is 
shown in Figure 12. The symbols used are listed in Table 3- 



Fleet Tankers or Barges 




Tanker or Pipeline 



Figure 12 
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Table 5 

Explanation of Symbols 

Day of year being considered 
Demand for product, day i 
Input of product, day i 

Transfer of product from working to storage tanks, 
day i 

Inventory in working tanks, day i 
Invent pry in reserve tanks, day i 

Maximum allowable inventory in reserve tanks, day i 

Cargo size of -tanker, day i 

th 

Cargo size of the j barge, day i 
Number of tankers, day i 
Number of barges, day i 
Expec.ted number of tankers 
Expected number of barges 
Average cargo of tankers 
Average cargo of barges 
Maximum allowable transfer rate 
Maximum input 
Planned input 

Upper decision limit for reducing input 

Upper decision limit for transfer to reserve storage 

Lower decision limit for transfer from reserve 
storage 

Lower decision limit for input increase 
Capacity of working tanks 
Capacity of reserve storage tanks 
Random variable 
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Analysis of the Problem 

The requirement or demand for product on day, was 

considered a function of the inventory in the tanks, planned 
input, and a random variable representing the variations in 
the arrival of tankers. 

The input to the system on day i is designated as P^ 

and is considered a function of the inventory in the working 

tanks and the product to be transferred to the reserve stor- 
age tanks. T i , the product transferred to reserve storage 
on day i, is related to the inventory in both the working 
and reserve storage tanks. 

A few of the functions developed by Dickson for P^, T^ 

and are listed here simply to indicate the nature of some 

of the restrictions utilized in the model. 

Daily Input, P i 

P. = P Daily input is equal to planned input if 

XT 

the following conditions exist: 

(1) < C c ; i.e. inventory in the working 

tanks must be between the limits beyond which it becomes 
♦ 

necessary to increase or decrease input. 

(2) Ip^ + Pp ” Tp < F; i.e. planned input will not 
overflow the tanks after transferring product to reserve 
storage. 

Daily Transfer, T i 



T = T 
x i x i max* 



Transfer to reserve storage at 
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maximum transfer rate when the following conditions exist: 



(l) I_,^ ^ C^; i.e. when the inventory in the shipping 

tank exceeds the upper control limit. 



(2) S* i ^ T max + T si ; i.e. pumping T max barrels into 
tanks containing T . barrels will not exceed the allowable 

S J_ 

inventory S'^. 

Daily Demand 

The functions defining demand involve the arrival rate 
of tankers and the size of the cargo lifted. • In Dickson's 
model the arrival fate of tankers was assumed to be a 
Poisson distribution after study of actual arrival records. 
A provision was included which was designed to consider the 
large tankers as scheduled within certain limits but which 
would require the arrivals of small tankers to follow the 
above mentioned Poisson distribution. 



Solution Method 

With the foregoing functions defined, Dickson indicates 
that inventories on each successive day can be obtained 
from material balances. Having this information and the 
delay time on tankers, operating costs are calculated. 

The values sought from the: model of the system were 
the optimal volume of shipping or working tankage and the 
control limits: C and C which are the upper limits on re- 

C 

ducing input and upper limits for transfer to storage; and 
C L and C R which are the lower decision limits on transfer 
to storage and increasing input respectively. 
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The method employed to obtain a solution was the Monte 
Carlo technique which was used to simulate the arrival of 
tankers and permitted the model to be played rapidly on a 
computer for a period equivalent to about 40 years. The 
model was thus used to estimate operating costs while var- 
ious values were tried for some of the variables subject 
to control. By averaging the operating costs over several 
years, a comparatively good estimate of operating costs 
was anticipated. 

Results and Limitations 

The results obtained by Dickson contained a sizeable 
error which was attributed primarily to the assumption of 
Poisson arrival rates for the tankers. The possibility of 
an extreme number of tankers arriving on a particular day 
was eliminated by placing limits on the maximum possible 
value in the Poisson distribution for tanker arrivals. 
However, there was no limit on the number of consecutive 
days without an arrival. This error, and probably others, 
contributed to the sizeable error in operating costs. 

In reviewing the study, Dickson found that because 
tankage is normally constructed in large increments of 
50,000 to 100,000 barrels per tank, the exact value of the 
shipping tankage was not required and the solution was suf- 
ficiently accurate to allow selection of the optimum tank- 
age. The values of the decision limits, for transferring 
product to or from storage and increasing or decreasing in- 
put, were found to be within the required range of accuracy. 
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Dickson indicates that in any subsequent model, limits 
should be placed on the number of arrivals on any day and 
the number of days with no arrivals. He also indicates 
that future studies should include some type of "anticipator 
mechanism" since tanker arrivals are usually known two or 
three days in advance in actual cases. 

It is considered by the present author that this method 
is worthy of study and it is hoped that the indicated modi- 
fications would prove adequate to the military problem. It 
is not likely that a general model could be developed to 
apply to all terminal situations. The problem of incorrectly 
assumed Poisson distribution of tanker arrivals can surely 
be corrected. Tanker arrivals are probably best described 
as intermediate between Poisson and scheduled. There are 

indications that tanker arrivals follow more closely the 

68 

distribution pattern of servically correlated arrivals 



CHAPTER IV 



APPLICATION OP OPTIMIZATION 
TECHNIQUES TO MILITARY PETROLEUM PURCHASING POLICIES 

As the nation's largest single purchaser and consumer 
of petroleum products, the Defense Department is vitally 
interested in the problem of procurement at minimum cost. 
Total petroleum purchases by the Defense Petroleum Supply 
Agency were slightly in excess of one billion dollars dur- 
ing the fiscal year 1961. Relying principally on competi- 
tive bidding to secure low prices, the Defense Petroleum 
Supply Agency consistently obtains prices well below those 
in commercial markets. 

The discussion and analysis of purchasing policy pre- 
sented in this chapter is intended to develop an avenue of 
approach to some of the basic factors involved in competi- 
tive bidding in general and petroleum procurement in parti- 
cular. It is hoped that an analysis of such factors will 
serve to indicate those areas where some of the optimiza- 
tion techniques presented in this thesis may be utilized. 

Competitive Bidding 

The objective of competitive bidding is to purchase a 
specified product at minimum cost. In the process of pre- 
paring a bid, a potential contractor must exercise his 
strategy in competition with that of his rivals. Although 
the bidder carries a cushion of profit in his bid, the bid 
price is, nevertheless, tied to his costs. 
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If because of some technological development all 
bidders are able to reduce their costs, it is reasonable 
to assume that bid prices would also be reduced as long 
as sufficient competition existed. Similarly, if some of 
the restrictions placed upon the bidder by the purchaser 
are altered, a corresponding change in bid prices may be 
expected. Some of the restrictions placed upon the bidder 
in the conditions set forth in the invitation for bids are: 
the type of contract, delivery requirements, specifications, 
quantity, etc. These do not represent all of the constraints 
which the bidder must consider, but they are among the major 
restrictions established by the purchaser. These conditions 
set forth in the invitation for bids, necessarily have a con- 
siderable and perhaps measurable effect on the ultimate 
price paid by the purchaser. 

Analysis of Constraints 

In seeking to minimize ultimate purchase prices, the 
objective of the buyer becomes one of choosing the constraints 
set forth in the invitation for bids in such a manner as to 
fulfill the necessary requirements and to permit the bidders 
to minimize their costs. 

While there are a number of conditions set forth in 
the invitation for bids which act as constraints, there are 
three which appear to be most important and which may be 
subject to some degree of control by the purchaser. 
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The three constraints to be discussed are: 

1. Material specifications 

2. Quantity 

3. Timing of Procurement 
Material Specifications 

It is recognized that specifications for items to be 
purchased are rarely subject to the control of the purchas- 
ing agency. Specifications, however, do represent a major 
part of the cost of any petroleum item. This becomes evi- 
dent when the price paid for each additional octane number 
is noted. The problem of optimizing specifications is not 
a simple one. Fuel specifications for an individual air- 
craft may be determined, but selecting an optimum fuel 
specification for a wide range of aircraft would be an 
extremely difficult linear programming problem. Any attempt 
to determine optimum specifications must consider not only 
the minimum standards for successful operation, but safety 
factors and the consequences of failure. Fortunately, most 
ships and aircraft are designed to fit existing fuels so 
the problem of optimum specifications may not be real. The 
analysis of optimum specifications is best left to persons 
qualified in this area. 

Quantity 

Military petroleum contracts, for the most part, are 
written as indefinite quantity contracts under which the 
government contracts for an estimated quantity but is only 
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obligated to order a token amount. Under such an arrange- 
ment, the government protects itself where fluctuations 
in requirements occur. The price increase, if any, which 
the government pays for the uncertainty resting with the 
contractor is small compared to the advantage gained. 

Usually the quantity ordered under the contract very nearly 
equals the amount initially estimated and probably no price 
differential exists. 

Since most military petroleum contracts provide for 
delivery on an as-required basis, and since the requirements 
of any given terminal or station are reasonably constant 
over an extended period, the total quantity contracted for 
is roughly proportional to the length of the contract period. 
The constraint placed on the bidders in terms of quantity 
can thus be said to be the length of the contract period. 

Optimum Length of Contract 

Numerous military petroleum contracts have a contract 
life of six months. This is believed to be relatively short. 
Before discussing possible methods of determining the optimum 
length of a contract period, let us consider the effects of 
a change in the length of the contract. 

If the contract period were lengthened, inventory costs 
for the government would remain essentially unchanged since 
we are considering that petroleum products are ordered only 
as required. Similarly, inventory costs for the refiner 
would probably not change appreciably since the volume per 
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unit time required by the government would remain substanti- 
ally unchanged. However, with longer term contracts, the 
refiner could conceivably reduce his costs in the areas of 
equipment utilization, production scheduling, and possibly 
in the cost of raw materials. If our assumption is correct, 
reduced costs to the refiner should result in lower prices 
to the government. 

Let us consider the position of a refiner bidding on 
one million barrels of 115/1^5 aviation gasoline. To make 
this product requires use of an alkylation unit which will 
be ulilized 'at, perhaps, 80$ of capacity if a contract is 
secured but at only ^> 0 % of capacity without a contract. 

On the basis of a one year contract it may be that the re- 
finer can operate the alkylation unit to advantage, while 
with only a six month contract his costs must be spread 
over the shorter period. To fulfill a short term contract, 
the refiner may even find it cheaper to buy alkylate from 
another source and transport it to his refinery. In either 
event, the cost to the refiner appears to increase under 
the short term contract. 

Not all of the factors to be considered indicate longer 
term contracts are in order. The fluctuations in the market 
in terms of demand and price may cause commercial prices to 
drop while the government still holds a contract at a higher 
price. Such situations cannot always be predicted. It is 
also true that bidders sometimes shave their prices under 
that of the successful bidder on the previous contract. Such 
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practises are not to be expected over an extended period, 
however, and the refiners cost must ultimately be the basis 
of his bid. 

The solution of such problems must await the develop- 
ment of models, collection of data and application of a 
great deal of theory. Linear programming models of refinery 
processes and the refining industry have been developed as 
has been shown. Such models contain assumptions which make 
the determination of an optimum length of contract inaccur- 
ate. Nonetheless, it is considered that an indication of 
the direction in which to move could be determined. 

Timing of Procurement 

Certain fuels such as heating oil are subject to seasonal 
demands which result in fluctuations in market prices. In 
such cases, the timing of a solicitation for bids may be 
important. If purchases can be made at times when demand 
and price are at their seasonal lows, some saving can normally 
be expected. No extended discussion of seasonal purchasing 
will be attempted here, however, it is considered pertinent 
to note that the optimal time to purchase one product is 
not likely to be optimal for another product. 

Use of Optimization Techniques 

All of the foregoing is best reduced to quantitative 
terms. The solution. of the problems presented would involve 
the use of several of the basic techniques of linear pro- 
gramming, inventory theory, game theory, and others not 
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discussed in this thesis. Existing linear programming 
models of refinery systems might be adapted for use in 
determining solutions. Extensive data would be required 
and the cost of such an effort would be substantial. 

When it is realized that a reduction in the average 
purchase price of military petroleum products by as little 
as one per cent would result in an annual cost reduction 
of about ten million dollars, the potential of such a pro- 
ject becomes apparent. 

The thoughts expressed in this chapter are not those 
of an authority in the field of military petroleum procure- 
ment, and there may be ramifications not known to this 
author. Much additional study is required before any such 
project could be initiated. 
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CHAPTER V 



SUMMARY AND CONCLUSION 

In this thesis an attempt has been made to survey 
some of the modern analytical and objective methods for 
optimizing decisions. The optimization techniques pre- 
sented are applicable to problems encountered in the 
operation of military petroleum systems. 

The specific techniques covered include: linear pro- 
gramming, game theory, and inventory theory. Certain as- 
pects of queueing theory and the Monte Carlo method were 
described briefly in connection with a problem involving 
optimization of tanker-terminal systems. 

Sample problems in the area of military petroleum 
management were developed and solved utilizing some of 
the methods presented. Problems known to be of current 
concern to military petroleum management, were discussed 
with the intention of developing approaches to possible 
solutions. 

Time was not available to develop computer programs 
for the methods presented. While some problems can be 
solved by simple arithmetic, it is considered that use of 
a digital computer is necessary to make significant use 
of the more sophisticated techniques such as linear pro- 
gramming . 
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Some of the problems considered represent extensions 
of existing methods which could be applied directly or 
which may be found in current use by military supply acti- 
vities. Others will require additional investigation be- 
fore the problem can be precisely formulated. 

The necessity for precise problem formulation should 
be emphasized for it is central to the philosophy of modern 
methods for the optimization of operational systems. One 
is dealing with problems which must be reduced to precise 
mathematical form, and a great amount of ingenuity and 
creativeness is necessary to translate a problem into the 
terms of a meaningful mathematical model. 

The task is difficult on two accounts. First, one 
must express ideas in mathematical terms which heretofore 
have not been regarded as quantitative. Second, the formu- 
lation must be in a form sufficiently tractable that a 
numerical solution is feasible. It is for this reason that 
a large part of the effort of this thesis was devoted to 
problem formulation and the development of analogies between 
systems which have been analyzed successfully and systems as 
yet unanalyzed within the military petroleum supply system. 

In conclusion, it is hoped that the material contained 
herein will be of value to those concerned with the manage- 
ment and analysis of the operation of the military petroleum 
supply system, and will serve as a basis for further study. 
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